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PREFACE 


The  analysis,  computer  programs,  and  results  discussed  herein  are  the 
result  of  analytical  research  performed  for: 


Air  Force  Geophysics  Laboratory 

Hanscom  Air  Force  Base,  Massachusetts  01731 


I.  INTRODUCTION 


A  large  portion  of  the  work  performed  in  1984  under  the  provisions  of 
contract  FI 9628-84-C-0071  consisted  of  developing  a  model  for  infrared  emission 
from  CO^  molecules  in  the  upper  atmosphere.  This  effort  included  the  develop¬ 
ment  and  refinement  of  computer  codes  having  the  purpose  of  1)  calculating  the 
rate  coefficients  for  absorption  of  solar  flux  in  the  upper  atmosphere;  2) 
calculating  the  steady-state  populations  of  vibrational  levels  of  CO^  as  a 
function  of  altitude  for  specified  conditions  by  equating  the  excitation  and 
de-excitation  rates  of  all  important  collisional  and  radiative  processes;  and 
3)  predicting  the  radiance  to  be  observed  by  hypothetical  sensors  in  limb-  and 
zenith-look  geometry.  These  codes  have  been  applied  to  a  special  problem, 
namely  determining  the  limb  radiance  in  the  2. 5-2. 9pm  wavelength  region  due  to 
emission  in  the  101  combination  bands  of  CO^  during  daylight  hours.  In  our 
first  three  quarterly  reports  we  have  outlined  the  principal  physical  results 
of  this  study;  in  Section  II  of  this  final  report  we  present  a  more  detailed 
version  of  our  methods  and  conclusions. 

The  absorption  and  radiance  codes,  although  applied  to  the  C0o  radiance 
problem,  were  developed  for  use  with  other  absorber  and  emitter  molecules  as 
well;  in  fact,  they  are  quite  general  line-by-line  codes. 

During  the  last  quarter  we  improved  the  Nl.TE  algorithm,  which  is  the  basic 
radiance  calculation,  by  making  it  more  efficient.  It  is  now  fast  enough,  for 
all  CC>2  bands,  to  run  in  iterative  applications  where  it  may  be  invoked  many 
times  in  a  single  job.  Some  of  these  improvements  are  discussed  in  Section  111. 

The  other  major  emphasis  of  the  fourth  quarter  was  on  the  A.F.G.L.  Auroral 
Model.  The  codes  comprising  this  model,  known  collectively  as  the  Auroral 
Atmospheric  Radiance  Code  (AARC) ,  were  extensively  revised  and  made  compatible 
with  each  other  so  that  they  could  be  used  in  a  single  interactive  application. 
In  Section  IV  we  describe  the  calculations  performed  by  each  component  of  AARC 
and  the  modifications  which  were  made. 
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The  other  tasks  carried  out  during  the  life  of  this  contract  were  described 
in  the  first  two  quarterly  reports.  Notable  among  these  is  program  NLTEA  for 
auroral  radiance,  which  will  be  incorporated  into  AARC.  This  task  was  summa¬ 
rized  in  the  second  quarterly. 


II.  METHODS  AND  RESULTS  OF  THE  C02  2.7pM  MODELLING  EFFORT 
1.  INTRODUCTION 

In  order  to  learn  more  about  emission  in  the  infrared  from  non-LTE  regions 
of  the  atmosphere,  and  in  order  to  provide  a  vehicle  for  developing  a  compre¬ 
hensive  line-by-line  model  for  CO^*  we  have  carried  out  a  detailed  set  of  cal¬ 
culations  of  C02  limb  radiance  in  the  2.7pm  wavelength  region  for  daylit 
atmospheres.  Under  the  assumed  conditions,  which  simulate  those  of  the  daylit 
scans  of  the  SPIRE  experiment  <l-3>,  the  principal  physical  processes  involved 
are  the  absorption  of  solar  flux  by  high-lying  ro-vibrational  states  of  C02 
molecules  at  wavelengths  of  2.7pm  or  less;  the  collisional  redistribution  of 
populations  among  the  high-lying  vibrational  levels;  spontaneous  emission  in 
the  101  combination  bands  and  in  other  bands;  and  exchange  of  vibrational 
quanta  with  N  .  In  this  section  we  describe  calculations  which,  taken  together, 
constitute  a  model  for  infrared  emission  from  C02  at  2.7pm.  We  also  discuss 
some  of  the  implications  of  our  results  and  also  of  the  comparison  we  make  with 
the  experimental  results.  We  do  not  discuss  the  codes  in  detail,  but  we  note 
that  considerable  attention  has  been  paid  to  making  them  general,  flexible,  and 
efficient  so  that  they  may  be  applied  to  situations  other  than  those  which  are 
discussed  below. 

An  earlier  study  <4>  has  found  that  the  daytime  emission  detected  by  the 
SPIRE  sensor  near  2.7pm  contains,  primarily,  contributions  from  the  Vj  and  v3 
fundamentals  of  water  vapor  and  the  101  combination  bands  of  CC^.  In  the  latter 
case  it  was  also  established  that  solar  pumping  of  the  201  and  301  levels 
followed  by  fluorescence  to  100  and  200,  respectively,  gives  a  more  significant 
component  of  the  C02  emission  than  the  resonant  101 — 000  transitions.  The 
significance  of  the  201-100  transitions,  and  other  conclusions  of  the  earlier 
work,  are  confirmed  in  the  present  effort.  Additional  contributions  from  other 
bands  are  evaluated  as  well.  Several  collisional  excitation  and  quenching 
mechanisms  which  were  not  previously  considered  have  been  included  here,  and 
their  effect  is  assessed.  The  inclusion  of  these  collisional  mechanisms  de¬ 
creases  the  observed  emission  by  only  about  30%.  However  it  markedly  changes 


the  spectral  distribution  of  the  observed  radiance.  We  also  compare  our  pre¬ 
dictions  with  the  SPIRE  dataset.  This  comparison,  made  in  conjunction  with 
certain  trial  modifications  of  our  basic  model,  allows  us  to  put  limits  on 
the  values  of  some  of  the  poorly-known  collisional  rate  constants  entering  our 
calculations . 


2 *  CENF-RAL  approach  to  calculation  of  vibrational-level  populations 

Before  the  2.7pm  emission  from  CC^  can  be  computed,  the  steady-state  popu¬ 
lations  of  all  the  radiating  levels  must  be  calculated.  Figure  1  is  an  energy 
level  diagram  for  the  principal  (626)  isotope  showing  the  vibrational  states 
which  are  necessary  for  this  calculation.  The  2.7pm  transitions,  including 
solar  pumping  from  the  ground  state,  are  indicated  by  the  solid  lines.  The 
important  solar  pumping  at  2.0  and  1.6pm  is  also  indicated.  Four  groups  of 
high-lying  levels,  including  all  the  levels  radiating  at  2.7pm,  are  identified. 
Within  each  group  the  vibrational  levels  are  strongly  coupled  by  collisional 
interactions.  The  principal  radiating  levels  are  lOOlx  (x=l,2)  in  Group  1; 
llllx  (x=l,2)  in  Group  2;  2001y  (y=l,2,3)  in  Group  3;  and  3001z  (z=l,2,3,4)  in 
Group  4.  The  notation  for  the  vibrational  levels  is  standard  AFGL  notation  <5>, 
as  is  the  isotopic  designation. 

Table  1  lists  the  26  bands  in  the  2.7pm  region  whose  contribution  to  the 
total  radiance  has  been  calculated  in  this  modelling  effort.  These  bands  are 
divided  into  eight  sets,  A  through  G,  to  facilitate  our  discussion. 

The  only  mechanisms  which  couple  states  within  any  single  group  are  col¬ 
lisional  interactions  which  cause  small  changes  in  the  internal  energy  of  the 
CO^  molecule  and  reorient  the  vibrational  motion.  These  result  from  col¬ 
lisions  with  air  molecules,  whose  internal  energy  states  remain  unchanged. 
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ENERGY 


Figure  1.  C02  energy  levels  relevant  for  the  2.7pm  radiance  calculations,  in¬ 

cluding  the  four  groups  of  high-lying  states  discussed  in  the  text. 
Solid  lines  show  2. 7pm  radiative  transitions.  Vertical  lines  show 
the  principal  solar  pumping  channels.  For  Group  2,  the  2 . 7. :m 
resonance  bands  are  the  major  channels.  Channels  for  quenching  by 
V-V  interactions  with  N?  are  not  shown. 
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Vib.  Transition 

A 

Band 

Isotope 

(cm  1 ) 

(yn) 

( sec  1 

10012-00001 

626 

3612.8 

2.77 

11.0 

10011-00001 

626 

3714.8 

2.69 

16.9 

20013-10002 

626 

3568.2 

2.80 

18 . 1 

20012-10002 

626 

3692.4 

2.71 

24.4 

20012-10001 

626 

3589.7 

2.79 

16.0 

20011-10001 

626 

3711 . 5 

2.69 

33.5 

11112-01101 

626 

3580.3 

2.79 

10.8 

11111-01101 

626 

3723.2 

2.69 

16 . 5 

12212-02201 

626 

3552.9 

2.81 

10.6 

12211-02201 

626 

3726.7 

2.68 

18 . 1 

30014-20003 

6  26 

3527.6 

2.83 

25.2 

30013-20003 

626 

3679.6 

2.72 

26.1 

30013-20002 

626 

3556.8 

2.81 

28.3 

30012-20002 

626 

3676.7 

2.72 

43.9 

30012-20001 

626 

3550.7 

2.82 

16.1 

30011-20001 

626 

3705.9 

2.70 

49.6 

22213-12202 

626 

3518.7 

2.84 

20.0 

22212-12202 

626 

3703.5 

2.70 

24 . 1 

22212-12201 

626 

3527.8 

2.83 

12.8 

22211-12201 

626 

3713.8 

2.69 

41.4 

10012-00001 

636 

3527.7 

2.83 

8 . 7 

10011-00001 

636 

3632.9 

2.75 

15.7 

10012-00001 

628 

3571 . 1 

2.80 

13.5 

10011-00001 

628 

3675.1 

2.72 

13.1 

10012-00001 

627 

3590.9 

2.78 

12.1 

10011-00001 

627 

3693.6 

2.71 

15.6 

Table  #1  C02  vibrational  bands  whose  contributions  to  the  total 

radiance  observed  near  2 . 7  um  have  been  calculated  in 
this  study.  Bands  are  grouped  into  sets  A  through  G 
to  facilitate  the  discussion. 
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The  individual  states  within  each  group  are  also  coupled  to  other  (lower 
lying)  C02  states  by  the  absorption  and  emission  of  radiation,  mostly  at  2.7pm 
and  4.3pm,  and  by  exchange  of  vibrational  quanta  with  N2  molecules  <6-8>. 

The  latter  process,  the  net  effect  of  which  is  usually  to  quench  C02  2.7pm 
emission,  involves  a  change  of  ±1  in  the  V3  quantum  number.  In  general,  then, 
radiative  transitions  near  4.3pm  are  in  parallel  with  collision-induced  V-V 
transitions  coupling  the  same  vibrational  levels.  The  populations  of  the 
lower-lying  states,  which  are  needed  to  calculate  excitation  rates,  are  taken 
to  be  static.  This  is  a  reasonable  approximation,  as  these  levels  are  more 
strongly  coupled  to  each  other  than  to  the  higher  lying  levels  in  which  we 
are  primarily  interested.  The  static  populations  that  we  used  were  obtained 
through  the  courtesy  of  T.  Degges  of  Visidyne  Corp.  from  a  simulation  of  the 
conditions  of  the  SPIRE  experiment  using  the  A.F.G.L.  High-Altitude  Infrared 
Radiance  Model  (IRR)<9>.  We  shall  see  that  the  lower-lying  excited  states 
are  important  for  our  calculations  only  at  altitudes  below  60  km.  In  this 
altitude  range  our  static  populations  are  practically  the  same  as  LTE  popu¬ 
lations. 

The  four  groups  of  states  are  treated  as  independent  of  each  other.  This 
is  a  good  approximation  because  the  possible  transitions  between  states  in 
different  groups  are  very  much  weaker  than  other  transitions  involving  them. 
That  is,  because  of  larger  energy  differences,  collisional  transfer  rates 
are  lower  between  groups  than  within  groups.  The  most  important  radiative 
transitions  between  groups,  v2  transitions,  turn  out  to  be  a  minor  pertur¬ 
bation  in  the  overall  picture. 

In  order  to  calculate  the  steady-state  populations  of  vibrational  levels 
within  a  group,  one  simply  equates  the  sum  of  excitation  rates  for  each  level 
due  to  all  important  mechanisms  with  the  sum  of  de-excitation  rates  for  that 
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Reaction 


Forward  Rate  Constant 


1.  CO2(10011)  + 

2.  CC>2  (02211)  + 

3.  CO2<10011)  + 


M  t  C02 (10012)  +  M 
M  Z  C02 (10012)  +  M 
M  Z  CO2(02211)  +  M 


k13  =  3x10 
k12  =  1.5x10 

k23  =  2kl  2 


4. 

C02 (10012) 

+ 

N2 

(0) 

5. 

C02 (02211) 

+ 

N2 

(0) 

6. 

C02  (10011) 

+ 

N2 

(0) 

C02 (10002) 

+ 

n2(1)‘> 

C02 (02201) 

+ 

n2(1) 

co2(ioooi) 

+ 

n2(D. 

k  =  5x10 
v 


—  13/  300' 


Table  #2  Reactions  resulting  in  non-radiative  transitions  involving 
the  three  states  of  Group  1.  The  forward  rate  constants 
are  also  listed. 
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level.  We  illustrate  this  by  considering  the  10012  state  alone,  and  then 
a  group  of  states  (10012,  02211,  10011)  which  were  treated  simultaneously. 


In  the  simplest  view,  10012  is  populated  by  solar  pumping  from  the  ground 
state,  quenched  by  collisions  with  ground-state  nitrogen  and  additionally  de¬ 
pleted  by  spontaneous  radiation  in  the  4.3pm  and  2.7pm  bands.  If  we  allow 
the  subscript  "1"  to  stand  for  10012,  and  "G"  for  the  ground  state  denoted 
by  00001,  the  rates  for  these  three  processes  are  J  n  ,  k  [N2(0)]n,,  and 
A^n^,  respectively.  The  solution  for  the  desired  number  density  n^  is  found 
by  equating  the  first  rate  with  the  sum  of  the  second  and  third. 

In  these  expressions,  J.  is  the  number  of  photons/sec  absorbed  in  the 

1  (j 

band  00001-10012  per  molecule  of  C02,  n^  is  the  total  number  density  of 
C02,  Aj  is  the  total  Einstein  coefficient  for  10012,  is  the  rate  constant 
for  quenching  by  nitrogen,  and  [ N2 (0) ]  is  the  ground-state  population  of 
nitrogen.  All  these  quantities  are  presumed  known.  A^  has  a  value  of  about 
455  photons/molecule-sec,  with  a  branching  ratio  of  about  40:1  in  favor  of 
the  4.3pm  transition.  One  can  see,  then,  that  of  all  the  solar  energy  pumped 
into  the  10012  level,  very  little  re-emerges  as  radiation  at  2.77pm  even  at 
altitudes  where  collisional  quenching  is  negligible.  This  is  the  principal 
justification  for  neglecting  the  pumping  due  to  2.7pm  emission  from  C02  mole¬ 
cules  in  the  atmosphere. 

This  simple  picture  is  complicated  only  very  slightly  by  the  recognition 
that,  with  a  sufficiently  large  population  in  the  low-lying  10002  level  (which 
is  symbolized  by  the  subscript  A),  excitation  of  10012  may  occur  because  of 
collisions  with  excited-state  nitrogen  [ N 2 ( 1 ) ]  and  also  because  of  solar  pump¬ 
ing  the  4.3pm  band.  That  is,  another  solar  pumping  coefficient  J  must  be 

1 A 

determined,  and  additional  excitation  rates  J.  n  and  k'[N2(l)]n  must  be 

1 A  I  V  A 

included  in  the  equation. 

Because  of  the  expected  strong  collisional  coupling  of  the  10012  level 
with  the  02211  and  10011  levels,  this  simple  picture  is  inadequate  to  de¬ 
termine  the  population  of  10012.  One  must  identify  the  important  collisional 


interactions  (see  Table  2),  put  in  collisional  rate  constants,  and  solve  siraul 
taneous  equations  for  the  populations  of  all  three  levels.  Before  writing 
these  equations,  it  is  desirable  to  comment  on  the  effect  of  doing  this. 

First  of  all,  the  discussion,  above,  of  the  10012  level  also  applies 
to  10011  to  the  extent  that  one  can  ignore  collisional  coupling.  The  02211 
level,  however,  is  different  because  its  radiative  transitions  to  and  from 
the  ground  level,  and  to  and  from  01101,  are  forbidden.  As  a  result,  what 
population  there  is  comes  from  absorption  of  solar  photons  around  4.3pm  by 
the  thermally  populated  state  (02201)  at  1335  cm  1 .  The  pumping  is  there¬ 
fore  much  weaker  than  that  for  the  other  two  levels.  The  mechanisms  for 
quenching  emission  at  2.7pm,  on  the  other  hand,  consist  of  V-V  interactions 
with  nitrogen  molecules  and  radiation  at  4.3pm,  and  are  comparable  in  strength 
to  those  of  the  other  two  levels.  The  result,  then,  of  coupling  02211  to 
10012  and  10011  is  to  provide  a  pathway  whose  net  effect  is  to  deplete  the 
populations  of  the  latter  two  states.  (At  least,  this  is  true  for  altitudes 
at  which  a  sufficient  number  of  collisions  occur  to  couple  these  states  ef¬ 
fectively  . ) 

The  second  notable  effect  of  coupling  the  levels  within  this  group  is 
to  deplete  the  higher-lying  levels  (e.g.,  10011)  in  favor  of  the  lower-lying 
ones.  This  causes  a  larger  fraction  of  the  emitted  radiation  to  appear  at 

2.77pm - and  a  smaller  fraction  at  2.69pm - than  would  otherwise  be  the  case, 

and  thus  affects  the  distribution  of  radiation  in  the  wavelength  region  under 
consideration.  In  general,  the  levels  split  by  Fermi-resonance - levels  dif¬ 
fering  only  in  the  fifth  index - are  probably  coupled  to  each  other  much  more 

strongly  than  to  the  other  nearby  levels  such  as  02211,  which  have  a  different 
symmetry.  (We  call  these  other  levels  "side  levels"  for  want  of  a  better 
term.)  They  are  therefore  likely  to  be  in  equilibrium  with  each  other  at 
quite  high  altitudes,  considerably  higher  than  the  altitudes  at  which  they 
are  effectively  decoupled  from  the  "side  levels". 

The  equations  which  need  to  be  solved  for  the  number  densities  of  the 
three  states  are  written  in  matrix  form,  below.  The  matrix  elements  are 
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coefficients  derived  from  the  rate  constants  describing  the  collisional  inter¬ 
actions,  V-V  interactions  with  N2»  and  from  the  radiative  transition  proba¬ 
bilities.  The  subscripts  (1,2,3)  refer  to  (10012,  02211,  10011),  respectively, 
and  the  subscripts  (G,A,B,C)  refer  to  (00001,  10002,  02201,  10001).  J  is 
the  sum  of  all  solar  pumping  rates  for  level  i  (i  =  1 ,  2 ,  or  3) .  is  the 

sum  of  the  Einstein  coefficients  for  level  i,  and  the  K. .  are  the  rate  coef- 

ij 

ficients  for  reactions  coupling  states  i  and  j,  including  those  for  which  i 
and  j  are  within  the  same  group  and  those  for  which  i  or  j  represents  a  lower- 
lying  state.  The  diagonal  elements  of  the  matrix  are  the  de-excitation  rate 
coefficients  and  all  the  other  terms,  including  the  inhomogeneous  ones,  rep¬ 
resent  excitation  mechanisms.  The  result  is 


VKA1+K21+K31 

-K12  -K13 

nl' 

JlnT+KlAnA 

K21 

A2+KB2+K12+K32  ~K23 

n2 

= 

J2nT+K2BnB 

'Si 

“K32  A3'tKC3+K13"fK23 

"3 

,3nT4K3CnC 

(1) 


VC  .  5  .  , 


The  rate  coefficients  for  quenching  v^.e.,  ^ 

to  exchange  of  quanta  with  nitrogen  are  given  by 


kai  =  VN2(0)]  (2a) 


K1A  =  kv|N*(1)]  (2b) 


where  k  is  listed  in  Table  2.  The  reverse  rate  constant  is 
v 

k'  =  k  exp[-C2(E  -E  -2330)/T]  (3) 

V  V  1  A 

where  C2  =  1.4388  is  the  second  radiation  constant  and  T  is  the  kinetic  temper¬ 
ature  . 

The  rate  coefficients  for  intramolecular  exchange  of  quanta  are 

K ,  „  =  k,  „ [M]  K„ ,  =  k,\[M]  (4a) 

(4b) 
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II 

K31 
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K32 

k23 
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(4c) 


Upper 

State 

Lower  o 

State  § 

t— * 

11101 

10001 

12201 

10002 

1101 

13301 

1 102 

0001 

2201 

14401 

0002 

2202 

10003 

10012 

J 

J,v 

GROUP  1  02211 

J,V 

iJOll 

J 

J,V 

11112 

J 

J.v 

GROUP  2  03311 

J.v 

11111 

J 

J,V 

20013 

1 

J  ! 

J 

J.V 

12212 

J  I 

J 

J.V 

20012 

J  i 

J  J 

J.V 

GROUP  3 

i 

04411 

J  9  V 

12211 

j 

J 

J.V 

20011 

j 

J 

J.V 

30014 

j 

22213 

J 

30013 

j 

GROUP  4  22212 

J  J 

30012 

j 

22211 

J 

30011 

j 

TABLE  3.  Excitation  mechanisms  considered  in  the  calculation  of  the  pop¬ 
ulations  of  the  vibrational  levels.  "J"  stands  for  solar  pumping  and  "V" 
stands  for  excitation  by  exchange  of  vibrational  quanta  with  N2.  The  im¬ 
portant  de-excitation  mechanisms  are  just  the  processes  which  are  comple¬ 
mentary  to  these,  except  that  the  seven  levels  in  Group  4  are  also  quenched 
by  N2  collisions  and  spontaneous  emission  at  2.7um  and  4.3am.  Collisions] 
interactions  within  each  group  are  not  listed. 
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where  [M]  is  the  number  density  of  air  molecules  and  the  forward  rate  constants 
are  given  in  Table  2.  Since  we  have  written  the  forward  rate  constants  with 
two  subscripts,  with  the  first  identifying  the  lower  state  and  the  s  cond 
the  upper  state,  the  reverse  rate  constants  are  given  by 


8i 

kT.  =  k..  —  expI-C2(E.-E.)/T] 


( 


where  the  g's  are  the  statistical  weights  and  the  E's  are  the  vibrational 
energies  in  wavenumbers. 

The  three  other  groups  of  states  are  treated  in  ways  which  are  analagous 
to  that  described  above.  Table  3  lists  the  mechanisms  by  which  the  various 
levels  in  each  group  are  coupled  to  lower-lying  levels.  In  addition,  each 
level  in  each  group  is  coupled  to  all  the  other  levels  in  its  own  group.  In 
both  cases,  collisional  processes  are  repr  .*sented  by  rate  constants  which 
are  identical  or  analagous  to  those  in  Table  2  and  rate  coefficients  like 
those  in  Eqs.  (2)  and  (4).  The  precise  energy  levels  can  be  found  in  Refer¬ 
ence  11,  as  can  the  band  strengths  which  were  used  to  calculate  the  Einstein 
coefficients  for  all  possible  spontaneous  transitions. 

In  the  next  two  sections  we  give  some  of  the  results  of  the  solar-pumping 
and  vibrational-population  calculations,  and  indicate  the  relative  importance 
of  some  of  the  factors  we  have  considered. 


3.  SOLAR  PUMPING  CALCULATIONS 

A  line-by-line  code,  SABS,  which  uses  data  collected  in  the  A.F.G.L. 
Atmospheric  Absorption  Line  Parameters  Compilation  <5,  12>  computes  the  ab¬ 
sorption  of  solar  flux  as  a  function  of  altitude.  We  used  a  standard  solar- 
flux  profile  <13>  and  a  solar  elevation  angle  of  12c  ,  which  closely  corresponds 
to  the  conditions  of  the  first  daylight  scan  of  the  SPIRE  experiment.  Given 
the  temperature  and  total  C02  profiles,  the  program  calculates  the  trans¬ 
mission  through  and  absorption  within  shells  one  kilometer  thick  starting 


at  the  top  of  the  atmosphere  and  proceeding  downward  to  an  altitude  of  40  km. 
The  program  assumes  all  atmospheric  parameters  to  be  constant  within  each 
shell,  uses  a  rotational  temperature  equal  to  the  kinetic  temperature,  and 
represents  the  shape  of  each  absorption  cross-section  with  a  Voigt  function. 

The  absorption  profiles  at  each  altitude  are  then  found  at  a  sufficient  number 
of  points  within  each  line  to  enable  integration  over  frequency.  Summing 
the  integrated  absorption  over  lines  gives  the  number  of  transitions  per  second 
per  molecule  of  C02  between  the  lower  vibrational  state,  i,  and  the  upper 
state,  j - that  is,  the  coefficient  J... 

ji 

Figures  2-5  give  this  absorption  as  a  function  of  altitude  for  those 
transitions  of  the  major  isotope  which  populate  states  in  the  four  groups 
under  consideration.  Some  of  the  curves  represent  negligible  contributions 
to  the  total  solar  excitation  of  their  respective  upper  states  over  most  or 
all  of  the  altitude  range  and  are  included  only  for  completeness.  Curves 
for  the  636  isotope  are  also  included  in  Figure  2. 

Figures  2,  4,  and  5  demonstrate  that  greatest  absorption  usually  occurs 
in  bands  originating  in  the  ground  state,  while  the  hot  bands  make  smaller 
contributions.  However,  it  is  important  to  note  that  hot  bands  do  not  always 
give  negligible  contributions.  They  also  illustrate  the  point  that  the  shape 
of  the  curves  for  the  hot  bands  reflects  the  assumed  lower-state  population 
more  than  anything  else.  In  Figure  4,  for  example,  the  curves  representing 
absorption  from  the  lOOOx  levels  have  maxima  and  minima  at  the  same  altitudes 
as  the  corresponding  vibrational-temperature  profiles,  and  mimic  them  in  other 
ways  as  well.  (We  reiterate  that  in  the  low-altitude  regions  where  these 
bands  are  important,  the  vibrational  temperatures  nearly  equal  the  kinetic 
temperature.)  These  generalizations  are  of  course  only  vr’id  where  the  in¬ 
cident  solar  radiation  has  not  been  bleached  out  at  higher  altitudes. 


The  absorption  which  occurs  along  the  paths  used  to  calculate  these  curves 
ranges  from  negligible  ("thin"  lines)  to  severe  ("thick"  lines)  to  extreme. 

A  convenient  figure  of  merit,  t  ,  is  the  optical  depth  (at  a  given  altitude) 
at  the  center  of  the  strongest  line  of  a  band,  from  which  one  can  crudely 
estimate  the  attenuation  by  comparing  its  value  with  unity. 
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Figure  2. 


ALTITUDE  (  km) 

Solar  flux  absorption  coefficients  for  bands  populating  the  states  in 
Group  1 ,  as  a  function  of  altitude.  Solid  lines  refer  to  bands  of  the 
626  isotope,  dashed  lines  to  those  of  the  636  isotopes.  Xs  mark  high- 
altitude  limits  for  absorption  in  the  2.7um  bands  of  the  627  and  628 
isotopes.  Hot-band  curves  are  sensitive  to  the  static  population 
assumed  for  the  lower  states,  but  these  are  nearly  LTE  for  altitudes 
at  which  important  contributions  result. 
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Figure  4. 


Solar  flux  absorption  coefficients  for  bands  populating  the  states  in 
(roup  3,  ns  a  function  of  altitude.  All  curves  refer  to  the  626  isotope 


ABSORPTION  (photons/molecule  -  sec) 


The  two  uppermost  curves  in  Figure  2,  giving  the  transitions  up  from  ground 
to  the  lOOlx  levels,  exhibit  the  full  range  of  possibilities.  At  the  highest 
altitudes  the  curves  are  nearly  flat  because  the  C02  density  is  too  low  to 
significantly  deplete  the  incident  ridiation.  Between  100  and  55  km  the  solar 
excitation  drops  by  two  orders  of  magnitude  because  of  the  severe  absorption 
taking  place  in  this  region.  The  optical  depth  parameter,  t  ,  at  these  two 
altitudes  is  approximately  .1  and  200,  respectively.  The  extreme  case  is 
illustrated  by  the  upturn  near  50  km.  This  occurs  because  of  absorption  in 
the  wings  of  the  Voigt  lines,  which  become  more  prominent  because  of  collision¬ 
broadening.  At  these  altitudes  the  flux  at  the  line-centers  is  completely 
absorbed  out  but  because  of  the  small  collision  width  at  greater  altitudes, 
that  in  the  wings  is  not  diminished  at  all.  Thus,  the  increasing  breadth 
of  the  absorption  lineshape  is  responsible  for  the  increase  in  absorption, 
which  now  occurs  far  from  the  line  centers.  This  behavior  of  course  does 
not  occur  when  a  Doppler  linehape  is  used  to  describe  the  cross-sections; 
rather,  the  raonotonic  decrease  with  decreasing  altitude  is  maintained  below 
60  km  for  all  reasonable  temperature  profiles. 

For  the  other  bands  originating  in  the  ground  state  (c.f..  Figures  2, 

4,  and  5)  the  onset  of  severe  attenuation  occurs  at  lower  altitudes,  reflecting 

weaker  transitions,  or,  in  the  case  of  minor  isotopes,  smaller  ground-state 

populations.  In  Figure  2  the  absorption  curves  for  the  636  isotope,  which 

comprises  about  1.1%  of  the  total  C02,  are  flat  down  to  about  70  km,  where 

the  lines  begin  to  become  thick.  Curves  for  the  628  and  627  isotopes  are 

similar,  with  comparable  attenuation  at  about  60  and  47  km,  respectively. 

In  all  these  cases,  10%  attenuation  corresponds  to  values  of  about  .12  for 

T  .  The  high-altitude  limits  of  the  626  and  636  curves  are  in  practically 
c 

the  same  ratio  as  the  isotopic  abundances;  for  the  nonsvmmet r ical  isotopes 
this  is  also  nearly  true  despite  the  presence  of  twice  as  many  lines  in  each 
band . 

The  differences  between  the  10012  and  10011  curves  (Figure  2)  and,  more 
notably,  between  the  20012  curve  (top.  Figure  4)  and  the  others  in  the  00001- 
2001v  family,  occur  because  of  different  transition  strengths.  Attenuation 


of  solar  photons  has  the  effect  of  equalizing  the  pumping  rates  at  lower  al¬ 
titudes,  but  at  higher  altitudes  these  differences  strongly  affect  the  vibra¬ 
tional  level  populations  and  thus  the  distribution  of  radiation  in  the  2. 7pm 
window. 

As  mentioned  earlier,  the  absorption  in  the  hot  bands  is  most  strongly 
determined  by  the  populations  of  the  lower  states.  In  the  absence  of  at¬ 
tenuation,  the  absorption  curves  have  the  same  altitude-dependence  as  the 
vibrational-temperature  profiles  used  to  calculate  these  populations.  These 
profiles  all  have  maxima  at  about  49  km  and  most  of  the  hot-band  curves  in 
Figures  2-5  also  peak  at  this  altitude.  Exceptions  occur  for  the  bands  whose 
upper  states  are  in  Group  1  (Figure  2)  and  Group  2  (Figure  3),  where  the  maxi¬ 
mum  absorption  lies  above  60  km  because  of  the  attenuation  of  incident  flux. 

The  optical  depth  parameter,  t  ,  at  49  km  is  approximately  12,  12,  and  5  for 
these  three  cases,  indicating  that  absorption  is  significant  across  much  of 
the  bands.  In  fact,  the  sharp  dropoff  in  absorption  between  60  and  40  km 
occurs  more  because  of  the  attenuation  in  the  solar  flux  than  because  of  the 
decreasing  population  of  available  absorbers. 

The  effect  of  changing  the  solar  elevation  angle,  and  hence  the  distance 
traversed  through  the  atmosphere  by  sunlight,  is  important  only  when  attenua¬ 
tion  of  incident  flux  is  significant.  The  four  daylit  scans  of  SPIRE  actually 
correspond  to  angles  between  12°  and  6°.  Within  this  range  the  traversal 
distance  (-sec  0)  varies  by  a  factor  of  two.  Thus,  although  the  high-altitude 
limits  are  the  same  for  all  angles,  the  absorption  occurring  at  altitudes 
below  which  attenuation  of  incident  flux  is  significant  is  reduced  (by  compari¬ 
son  to  the  curves  given)  as  the  elevation  angle  decreases.  This  is  particularly 
important  for  the  00001 -lOOlx  bands  below  100  km,  and  for  the  00001-2001 v 
bands  (Figure  4)  below  75  km. 
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for  the  populations  of  the  levels  in  the  four  chosen  groups.  We  present  the 
results  as  the  ratio  of  the  population  in  each  vibrational  level  to  ground- 
state  C02  because  this  permits  the  easiest  graphical  visualization  of  the  ef¬ 
fect  of  subtle  but  important  changes.  We  also  present  curves  giving,  for  the 
most  important  vibrational  levels,  the  excitation  and  de-excitation  rates  (in 
units  of  transitions/sec-cm3)  as  a  function  of  altitude.  In  particular,  we 
break  out  the  separate  effects  of  radiative  excitation  and  de-excitation, 
and  also  the  net  effects  of  the  N2  vibrational-quantum  exchange,  the  col- 
lisional  interactions  connecting  the  Fermi-resonance-split  states,  and  the 
collisional  interactions  connecting  the  principal  radiating  levels  with  the 
side  levels. 

We  have  already  discussed,  qualitatively,  the  excitation  and  de-excitation 
mechanisms  responsible  for  the  populations  of  states  in  Group  1.  In  Figure  6 
we  give  the  population  ratios,  and  in  Figures  7  and  8  the  rates  for  10012 
and  10011,  respectively.  As  expected,  solar  pumping  is  dominant  for  popula¬ 
ting  these  two  levels.  Over  a  narrow  range  of  altitudes  near  50  km  the  net 
effect  of  N2  V-V  interactions  is  to  excite  these  states.  In  fact  the  prominent 
peaks  at  50  km  in  Figure  6  result  from  a  combination  of  this  excitation  by 
N2  and  the  solar  pumping  in  the  wings  of  the  Voigt  lines  which  was  discussed 
in  the  previous  section.  Between  54  and  80  km,  however,  N2  V-V  interactions 
constitute  a  strong  quenching  mechanism  which  is  largely  responsible  for  lim¬ 
iting  the  radiation  at  2.7ym  and  4.3um.  Only  above  75  km  does  the  spontaneous 
emission  become  the  most  important  mechanism  for  depopulating  these  levels. 

The  effect  on  the  two  lOOlx  levels  of  including  the  collisional  interac¬ 
tions  with  the  side  level  02211  is,  as  mentioned  earlier,  to  deplete  both. 

The  depletion  rate  can  be  seen  from  Figures  7  and  8  to  be  modest  but  not  neg¬ 
ligible.  Coupling  the  Fermi-resonance-split  levels,  on  the  other  hand,  is 
important  even  up  to  100  km,  where  its  net  effect  is  30%  of  the  solar  pumping 
rate . 

The  three  levels  comprising  Group  2 - 1  1  1  1  1  ,  03311  ,  1  1  1  12 differ  from 

those  of  Group  1  in  that  all  the  radiative  transitions  to  ground  are  very 
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Figure  6. 
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Figure  7.  Excitation  rates  (solid  lines)  and  de-excitation  rates  (dashed  lines) 
due  to  mechanisms  coupling  the  10012  level  with  other  levels.  "J" 
represents  solar  pumping  from  lower  states  (primarily  ground),  "E" 
represents  total  spontaneous  emission,  "VV"  represents  N;  collisions, 
and  10011  .and  "0221  1  "  represent  collisional  interactions  with  those 
two  level::. 


weak.  Nevertheless  there  is  strong  radiative  coupling  to  the  low-lying  01101 
state,  which  is  populated  strongly  enough  by  various  collisional  and  radiative 
processes  that  solar  pumping  to  the  111  lx  levels  is  very  important.  There 
is  also  radiative  and  N2-collision  V-V  coupling  to  the  11101  and  11102  levels, 
which  are  not  so  highly  populated.  The  net  effect  of  both  of  the  latter  is 
to  quench  the  higher  states,  but  in  the  altitude  range  below  55  km  it  is  im¬ 
portant  not  to  discard  the  excitation  component  of  these  processes  and  retain 
only  the  obvious  de-excitations.  Above  55  km,  however,  these  omissions  would 
be  of  little  consequence. 

Reference  to  Figures  9-11  shows  that,  apart  from  the  collisional  exchanges 
among  themselves,  11112  and  11111  are  populated  almost  entirely  by  solar  pumping 
Depletion  is  primarily  due  to  V-V  interactions  with  N2  below  74  km  and  spontan¬ 
eous  emission  above  this  altitude.  The  collisional  interactions  rearrange  the 
populations  within  the  group,  as  before.  The  collisional  coupling  to  the 
03311  side  level,  which  is  quenched  by  the  same  processes  but  lacks  an  im¬ 
portant  radiative  excitation  channel,  tends  to  deplete  the  other  two  levels. 

The  third  group  of  states,  including  the  2001y  levels,  is  simpler  to 
discuss  because  the  only  excitation  mechanism  of  any  importance  is  solar  pumping 
(The  net  effect  of  V-V  interactions  with  N2  is  so  strongly  to  quench  these 
states  that  it  makes  no  difference  if  this  is  disregarded  as  an  excitation 
mechanism.)  The  bulk  of  the  pumping  comes  from  ground  to  the  2001y  levels, 
but  near  50  km  there  are,  as  reference  to  Figure  4  shows,  important  contri¬ 
butions  due  to  absorption  near  2.7pm  both  to  these  main  levels  and  also  to 
the  12211  and  12212  side  levels.  Nevertheless  for  altitudes  above  60  km  the 
populations  of  the  six  states  can  be  explained  almost  completely  by  absorption 
of  2.0pm  solar  flux,  redistribution  by  collisional  interactions,  and  quench¬ 
ing  by  N2  at  lower  altitudes  and  by  radiation  at  higher  ones.  74  km  is  again 
the  altitude  at  which  these  quenching  mechanisms  are  of  equal  importance. 

As  is  typical,  the  side  levels  1221x  are  mostly  populated  by  energy  transfer 
from  the  2001y  levels.  Since  these  also  have  the  possibility  of  radiating 
at  2.7pm,  however,  including  them  does  not  greatly  diminish  the  total  emission 
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Ratios  of  the  populations  of  the  six  vibrational  levels  in  Groun 
to  that  of  the  y, round  state  for  the  626  isotope. 


29 


ifu  yu.A. 


& 


Wl' 

$ 

n'! 


w> 

c  S 
o  10° 


Figure  1 3. 


ALTITUDE  (km) 

Excitation  and  de-excitation  rates  for  the  20013  vibrational  level. 
Solid  lines  represent  net  excitation;  dashed  ones,  net  de-excitation 
"•!'  refers  to  solar  pumping,  "E"  to  spontaneous  emission,  "W"  to 
collisions  with  N -  molecules,  "F"  to  collisional  interactions  with 
other  Fermi- resonance  split  states  in  Group  3  (20011,20012),  and 
"S"  to  collisional  interactions  with  the  three  remaining  levels  in 
Group  3. 
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calculated  in  the  end,  even  if  there  is  some  redistribution  in  frequency. 
This  is,  however,  not  true  of  the  other  side  level,  04411. 


Figures  12-15  give  the  quantitative  results. 

Most  of  the  comments  used  to  describe  Group  3  also  apply  to  Group  4. 

In  this  case,  however,  the  solar  pumping  is  much  weaker  and  the  end  result 
is  populations  which  are  not  large  enough  to  contribute  appreciably  to  the 
total  observed  radiance.  We  plot  the  populations  of  these  levels  in  Figure  16 
for  the  sake  of  completeness. 


5 .  RADIANCE  CALCULATIONS 

The  radiance  calculations  were  performed  using  the  line-by-line  code 
NLTE  <14>  in  limb-look  geometry.  Most  of  the  procedures  are  similar  to  those 
mentioned  in  connection  with  program  SABS,  except  that  it  is  now  necessary 
to  account  for  continuous  emission  along  the  whole  1 ine-of-si ght  path. 

The  most  important  result  of  this  study  is  presented  in  Figure  17,  which 
gives  the  total  integrated  radiance  in  the  2.?um  region  as  a  function  of  the 
tangent  height  of  the  1 ine-of-sight  path,  and  also  breaks  it  down  into  contri¬ 
butions  from  the  sets  of  bands  identified  in  Table  1.  It  is  clear  from  these 
curves  that  the  dominant  contribution  at  high  altitudes  comes  from  the  reso¬ 
nant  transitions,  set  A,  while  at  low  altitudes  the  fluorescent  bands  of  set 
B,  bands  involving  transitions  2001v-1000x,  are  the  most  important  component. 
(Their  contributions  are  equal  for  the  path  with  a  tangent-height  of  approxi¬ 
mated  77  km.)  Next  in  importance  are  the  bands  of  set  C,  and  the  resonant  bands 
of  the  minor  isotopes,  summed  (set  G) .  Set  D  gives  10%  of  the  total  for  the 
60  km  path  but  drops  rapidly  in  importance  with  increasing  altitude,  as  do 
all  bands  whose  upper  states  are  populated  primarily  as  a  result  of  collisions. 
The  two  lowest  curves,  resulting  from  states  in  Group  4,  contribute  at  most 
27  of  the  total  and  need  not  be  considered.  In  addition,  the  2001y-1000x 
bands  of  the  minor  isotopes,  not  shown,  fall  in  the  range  between  these  last 
two  curves. 
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Figure  17 


75  80  85  y 

TANGENT  HEIGHT  (km) 

Band  radiance  in  the  2.6-2. 9um  range,  as  a  function  of  tangent  height. 
Contributions  from  the  different  sets  of  bands  identified  in  Table  1  are 
(7)  Set  A,  x  Set  B,  A  Set  C,  #  Set  E,  ^  Set  F,  o  Set  G. 


The  radiance  contributed  by  the  major- isotope  resonant  bands,  set  A, 
changes  very  little  over  the  range  of  line-of-sight  paths  because  self- 
absorpt ion  is  important  even  for  the  highest  paths.  On  the  100  km  path,  for 
which  these  bands  comprise  90%  of  the  total  radiance,  the  optical  depth  param¬ 
eters,  1  ,  are  about  1.3  and  2.0  for  10012-00001  and  10011-00001,  respectively, 
and  the  actual  contributions  are  greater  than  for  the  60  km  path.  All  the 
other  bands  have  only  thin  lines  on  the  100  km  path,  and  their  contributions 
simply  reflect  the  upper-state  populations,  which  of  course  decrease  rapidly 
with  increasing  altitude.  In  fact,  the  only  other  bands  for  which  self¬ 
absorption  plays  an  important  role  at  any  altitude  above  60  km  are  the  minor- 
isotope  resonance  bands,  set  G,  and  the  major-isotope  llllx-01101  hot  bands 
of  set  C,  whose  lower  state  is  only  667  crn-1  above  ground.  At  60  km  the 
largest  value  for  any  band  in  set  B  is  about  .5,  which  implies  a  little 
attenuation  along  that  path. 

One  can  see  from  Figure  17  that  between  60  and  75  km  it  in  important 
to  consider  at  least  five  contributions  to  the  total  radiance.  No  single 
family  comprises  more  than  45%  of  the  total  in  this  altitude  range.  On  the 
60  km  path  the  contributions  from  sets  B,  A,  C,  G,  and  D  are  45%,  16%,  13%, 

13%  and  11%  respectively;  for  the  75  km  path  they  are  37%,  28%,  15%,  14%  and 
5%. 

Figures  18  and  19  give  the  spectral  distribution  of  radiation  between 
~.6um  and  2 . 9um.  In  order  to  obtain  this,  the  NLTE  output  for  each  line  was 
convolved  with  an  instrumental  scanning  function  appropriate  to  the  SPIRE 
experiment.  The  lines  were  approximated  by  delta  functions  with  amplitudes 
equal  to  the  integrated  radiance.  A  triangular  scanning  function  with  a  FWHM 
of  52  cm-1,  or  .038um  <2-  was  used  to  simulate  the  effect  of  this  instrument. 

Figure  18  gives  the  total  predicted  spectral  radiance  from  all  the  C02 
bands  for  paths  with  tangent  heights  between  60  and  100  km.  One  notable 
feature  is  that  the  relative  size  of  the  two  peaks  changes  substantially  with 
tangent  height.  This  occurs  because  of  the  diminishing  importance,  with  in¬ 
creasing  altitude,  of  collisional  excitation  and  de-excitation  mechanisms 
as  compared  with  radiative  mechanisms.  That  is,  collisions  have  the  effect 
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Spectral  distribution  of  litnb  radiance,  from  the  basic  model, 
refer  to  tangent  heights  between  60  and  100  km. 
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of  diminishing  the  population  of  higher-lying  levels  within  each  group  and 
enhancing  that  of  the  lower  ones,  and  the  net  result  for  lower  altitudes  is 
more  radiation  at  longer  wavelengths.  On  the  other  hand,  at  100  km  the  radiative 
excitation  and  de-excitation  mechanisms  stand  alone  and  (because  the  pumping 
to  the  10011  level  is  more  likely  than  that  to  10012,  as  discussed  above) 
the  shorter-wavelength  peak  is  larger  for  paths  that  traverse  regions  near 
and  above  this  altitude. 

Another  point  worth  noting  is  the  shifting  of  the  1 onger-wavelength  peak 
from  about  2.80pm  to  about  2.77pm  as  the  tangent-height  increases  from  60 
to  100  km.  This  is  a  result  of  the  decreasing  importance  of  two  of  the  fluo¬ 
rescent  bands  in  set  B,  which  peak  at  slightly  longer  wavelengths,  when  com¬ 
pared  to  the  resonant  band  at  2.77pm.  The  shorter-wavelength  peaks  are  more 
nearly  coincident. 

Figure  19  gives  a  breakdown  of  the  total  spectral  radiance  seen  on  the 
75  km  path  by  sets  of  bands,  including  only  the  five  most  prominent  contri¬ 
butions.  It  is  interesting  to  note  that  the  long-wavelength  tail  of  the 
distribution  (\>2.85pm)  is  primarily  determined  by  the  minor-isotope  emission. 

6.  MODIFICATIONS  TO  THE  BASIC  MODEL 

In  order  to  assess  quantitatively  the  effects  of  including  the  intra¬ 
group  collisional  interactions,  we  performed  the  vibrational-population  and 
radiance  calculations  using  alternative  values  for  the  rate  constants. 

The  first  modification  (Modi)  consisted  of  setting  those  rate  constants 

which  couple  the  side  levels  to  the  principal  radiating  levels - e.g.,  k^9 

in  Table  2 - equal  to  zero.  This  has  the  effect  of  eliminating  the  influence 

of  these  levels  (02211  in  Group  1;  03311  in  Group  2;  12211,  12212,  and  04411 
in  Group  3;  22211-22213  in  Group  4)  on  the  vibrational  population  calculations. 

In  modification  number  2  (Mod2)  the  rate  constants  coupling  the  Kerin i- 
resonance  split  states - e.g.,  k^  in  Table  2 - as  well  as  those  coupling 


them  to  the  side  levels  were  set  equal  to  zero,  thereby  making  all  levels 
of  each  group  independent  of  all  others.  This  results  in  calculations  which 
are  most  like  those  previously  reported  <4>. 


Modification  number  3  (Mod3)  is  the  same  as  the  basic  model,  except  that 
the  rate  constants  connecting  the  Fermi-resonance  split  levels  were  increased 
by  a  factor  of  100.  These  levels  are  thereby  coupled  more  strongly  to  each 
other. 

Figures  20  and  21  summarize  the  overall  effect  of  making  these  alterations 
by  plotting,  respectively,  the  total  radiance  as  a  function  of  tangent  height 
and  the  2.7um  spectral  radiance  for  the  75  km  tangent-height  path,  each  cal¬ 
culated  under  the  assumptions  of  the  basic  model  and  those  of  Modi  and  Mod2. 
Mod3  produces  radiances  which  hardly  differ  from  those  of  the  basic  model. 

Two  things  become  clear  from  comparing  the  results  of  the  basic  model 
to  those  of  Modi  and  Mod2.  First  of  all,  the  rate  at  which  collisions  re¬ 
distribute  the  populations  of  the  Fermi-resonance  split  levels  among  them¬ 
selves  has  a  great  influence  not  only  on  the  spectral  distribution,  but  also 
on  the  total  emission.  Secondly,  the  decoupling  of  the  side  levels  has  the 
anticipated  effect  of  increasing  the  emission  from  the  principal  radiating 
levels,  but  it  is  bv  an  amount  which  is  considerably  less  than  that  re¬ 
sulting  from  decoupling  all  the  levels.  A  more  detailed  examination  of  the 
results  reveals  that  the  principal  differences  between  the  Mod2  curves  and 
the  others  comes  from  much  greater  populations  of,  and  hence  emission  from, 
certain  states  of  Group  3,  nrimarily  20012. 

The  most  general  conclusion  that  can  be  drawn  is  that  the  observed 
radiance  depends  very  strongly  on  the  collisional  interactions  among  the 
high-lying  states,  and  that  it  is  very  important  to  know  the  rate  constants 
for  the  reactions  coupling  these  states. 


To  see  how  the  changes  depicted  in  Figures  20  and  21  arise,  we  examine 
Modi  and  Mod2  in  more  detail.  Modi  decouples  the  side  levels  completely  and 
removes  a  quenching  mechanism  for  the  more  important  radiating  states.  In 
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Figure  6  the  dashed  line  gives  the  populations  of  Group  1  states  under  these 
assumptions.  It  shows  that  above  the  altitude  (approximately  55  km)  where 
N2  excitation  is  unimportant  the  02211  state  is  not  highly  populated.  The 
lOOlx  populations  are  augmented  by  as  much  as  20%  between  60  and  80  km,  but 
the  limb  radiance  resulting  from  them  is  increased  by  a  smaller  percent  be¬ 
cause  of  self-absorption.  In  Group  2,  the  111  lx  populations  are  augmented 
by  1j%  or  less. 

Modi  strongly  affects  the  populations  of  the  2001y  states  of  Group  3, 
with  increases  of  80%  at  some  altitudes  and  appreciable  effects  all  the  way 
from  40  to  90  km.  This  more  dramatic  response  to  the  decoupling  of  the  side 
levels  reflects  the  lower  excitation  rates  and  the  greater  number  of  quenching 
channels  (compared  to  the  lOOlx  levels)  for  these  states.  The  change  in  the 
total  limb  radiance  from  the  basic  model  is  mostly  due  to  these  augmented 
populations,  and  occurs  despite  the  near-complete  elimination  of  the  con¬ 
tribution  from  the  bands  of  set  D  originating  in  the  1221x  side  levels. 

The  Modi  spectral  distribution  is  not  significantly  different  from  that 
of  the  basic  model. 

The  second  modification - namely,  decoupling  all  the  high-lying  states 

from  each  other - has  very  important  implications.  In  Modi  the  Fermi-resonance 

split  states  are  still  strongly  coupled,  so  the  additional  limitation  imposed 
by  Mod2  results  only  in  the  redistribution  of  the  populations  of  the  ex¬ 
cited  states,  compared  to  Modi.  One  would  thus  expect  little  change  in  the 
total  emission  from  Modi  to  Mod2.  This  neglects,  however,  the  important 

fact  that  the  states  which  are  most  strongly  coupled  to  ground - and  which 

therefore  are  more  strongly  pumped - may  radiate  more  (or  less)  strongly  in 

their  respective  2.7ym  bands  than  other  states  in  the  same  group.  The  ab¬ 
sence  of  the  collisional  redistribution  will  therefore  result  in  more  (or 
less)  total  emission  from  these  states.  This  is  a  pertinent  consideration 
for  the  lOOlx  levels  and  the  2001v  levels.  In  fact,  the  Einstein  A  coef¬ 
ficient  for  the  transition  to  ground  of  the  10011  level  (16.9  sec-*)  is 
about  50%  greater  than  that  of  the  10012  level  (11.0  sec-*),  so  (as  can  be 
seen  from  Figure  2  above  90  km)  the  solar  pumping  of  this  level  is  greater 


and  also  the  probability  of  radiation  at  2.7yro  from  this  level  is  greater. 

Since  this  level  is  actually  depleted  in  favor  of  10012  by  collisional  inter¬ 
actions,  the  proper  result  is  less  emission  than  is  found  when  (as  in  Mod2) 
these  levels  are  assumed  to  be  independent. 

The  same  result  holds  for  the  2001y  family  of  Group  2.  One  can  see  from 
Figure  4  that  for  altitudes  at  which  attenuation  of  sunlight  is  unimportant 
the  20012  level  is  much  more  strongly  pumped  than  the  other  two.  It  happens 
that  the  total  Einstein  coefficient  for  emission  at  2.7pm  from  20012  is  greater 
than  that  of  the  lower-lying  20013  (18.1  sec-*)  and  higher-lying  20011  (33.5 
sec-l).  The  net  effect  of  decoupling  these  levels,  then,  is  to  keep  the  more 
weakly-radiating  lowest  level  from  acquiring  the  large  population  collecting 
in  the  middle  and  upper  levels.  Mod2  therefore  predicts  a  considerably  greater 
total  emission  near  2.7pm  than  the  basic  model,  as  shown  in  Figure  20.  This 
result  calls  attention  to  the  need  to  carefully  consider  the  collisional  inter¬ 
actions. 

Figure  21  shows  that  the  changes  which  constitute  Mod2  give  a  strikingly 
different  spectral  distribution  from  that  of  the  basic  model.  This  occurs 
because  of  the  greater  weighting  of  Mod2  accords  the  higher-lying  states, 
which  radiate  more  prominently  in  shorter-wavelength  bands. 

Mod3  is,  in  some  sense,  the  opposite  of  Modi  and  Mod2  because  it  couples 
the  principal  radiating  states  more,  rather  than  less,  strongly.  The  quanti¬ 
tative  result  is  that  the  populations  are  weighted  slightly  more  heavily  toward 
the  lower  vibrational  levels,  as  one  would  expect.  Below  85  km  the  population 
differences  are  too  small  to  be  seen  on  the  logarithmic  scale  of  Figures  2-5; 
above  85  km  the  redistribution  is  more  obvious.  By  the  converse  of  the  logic 
employed  in  the  discussion  of  Mod2  in  the  last  two  paragraphs,  one  expects 
this  to  reduce  the  observed  radiation.  This  turns  out  to  be  the  case,  but 
since  emission  from  Group  3  is  relatively  unimportant  above  85  km  the  limb 
radiance  is  not  affected  dramatically. 


7.  COMPARISON  WITH  EXPERIMENTAL  RESULTS 


The  2.7ym  emission  seen  by  SPIRE  <4-6>  and  other  experimental  probes 
has  components  due  to  molecules  other  than  C02.  The  most  important  of  these 
is  H20,  which  radiates  in  four  bands  between  2.68ym  and  2.75pm.  The  possi¬ 
bility  of  contributions  from  the  OH  fundamental  on  the  long-wavelength  side 
cannot  be  ruled  out,  either. 

Modelling  the  water  contribution  is  made  difficult  by  the  uncertainty 
in,  and  variability  of,  the  mixing  ratio  for  H20.  Degges  <15>  has  determined 
that  the  SPIRE  data  at  6.3ym  and  2.7um  for  paths  with  tangent  heights  up  to 
75  km  are  consistent  with  a  mixing  ratio  of  1.5  parts  per  million  by  volume, 
although  higher  values  are  not  ruled  out.  He  has  calculated  the  H20  contri¬ 
bution  in  the  2.5-2.9ym  region  on  this  basis.  In  Figure  22  we  compare  the 
results  of  our  basic  model  for  C02 ,  Degges'  results  for  H20,  and  the  total 
of  these  quantities  with  the  SPIRE  dataset  for  tangent  heights  between  50 
and  90  km. 

Figure  22  shows  that  below  70  km  there  is  good  agreement  between  the 
total  predicted  limb  radiance  and  the  experimental  values,  which  do  not  vary 
strongly  with  altitude.  Above  75  km  the  emission  falls  off  rapidly  into  the 
range  of  detector  noise.  In  fact,  the  data  are  very  noisy  at  these  altitudes, 
and  it  is  only  possible  to  say  that  the  actual  emission  must  be  less  than 
1.6xl0~®  watt/cm  -ster,  the  noise  level.  Since  the  C02  prediction  is  less 
than  this,  no  conclusion  can  be  drawn  regarding  it.  The  H20  prediction  of 
limb  radiance  is  significantly  above  this  for  tangent  heights  between  75  and 
85  km,  however,  a  fact  that  suggests  that  the  correct  water  vapor  mixing  ratio 
over  most  of  the  range  above  75  km  must  be  less  than  1.5  ppm. 

Figure  23  presents  spectral  data  in  the  2. 5-3. 5pm  region,  with  the  back¬ 
ground  subtracted  out.  These  data  correspond  to  a  tangent  height  of  73  km. 

We  also  show  the  H20  component  <15>  and  the  combined  C02  and  H20  prediction. 
The  latter  provides  a  respectable  representation  of  the  data.  It  somewhat 
overestimates  the  experiment  at  2.70um,  in  the  shorter  wavelength  C02  peak 
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(c.f.  Figure  18);  and  it  underestimates  it  on  the 
the  C02  spectrum  and  beyond.  This  fit  illustrates 
tance  of  including  the  collisional  interactions  in 
if  the  results  of  the  basic  model  were  replaced  by 
Figure  21)  both  the  overestimate  of  the  experiment 
estimate  beyond  2.80vim  would  be  exaggerated  rather 


long-wavelength  side  of 
quite  clearly  the  impor- 
the  C02  model.  That  is, 
the  Mod2  results  (c.f. 
at  2.70ym  and  the  under- 
than  diminished. 


It  should  be  pointed  out  that  the  water  calculation  is  a  band  model  cal¬ 
culation,  and  mav  not  represent  the  spectral  distribution  as  precisely  as 
appears  necessary  to  resolve  the  discrepancies  between  the  combined  model 
and  the  experiment  in  the  2.7-2.8um  region. 


8.  DISCUSSION 

It  is  clear  from  the  discussion  of  the  last  two  sections  that  the  inter¬ 
actions  coup! ing  the  states  within  each  of  the  groups  of  high-lying  states 
plav  a  major,  albeit  not  dominant,  role  in  determining  the  C02  contribution 
to  the  2.7um  radiance.  The  values  of  the  rate  constants  used  in  the  basic 
model,  given  in  Table  2,  wore  chosen  to  be  the  same  as  those  coupling  the 
10001-02201-10002  states  in  IRK  <9'.  In  fact,  is  the  lower  limit  pro¬ 

posed  by  Taylor  <10  -  f or  the  rate  constant  characterizing 

C02(  10001)  +  M  C02(  10002)  +  M 


Tlu  se  values  are  not  known  with  great  accuracy,  and  we  also  do  not  know  exactly 
what  differences  can  be  expected  between  the  constants  coupling  states  of 
this  low-lying  group  and  those  characterizing  our  four  high-lying  groups. 

The  values  used  in  our  basic  model  thus  represent  educated  guesses  of  how 
strong  the  interactions  are.  The  purpose  of  making  the  changes  described 
as  Modl-M(’d3  was  to  crudely  evaluate  the  effect  of  these  guesses  on  the 
observed  radiance. 


The  effect  of  reducing  kj  ,  can  be  seen,  from  Mod2,  to  be  to  increase 
the  2.7am  emission.  If  the  H20  contribution  predicted  bv  Degges  <15-'  is 
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correct,  the  overall  fit  in  the  50-70  km  region  precludes  the  possibility 
of  a  larger  C02  component.  Thus  the  value  we  have  used  for  k^  could  be  a 
lower  limit,  in  accord  with  Taylor  <10>.  We  have  not  shown  that  some  lower 
values  are  inconsistent  with  the  data,  but  since  k^  determines  the  largest 
de-excitation  rate  for  the  very  important  20012  and  20011  states  (c.f. 

Figures  14  and  15)  it  is  likely  that  modest  reductions  in  its  value  will 
lead  to  larger  populations,  and  therefore  greater  emission. 

Because  the  Fermi-resonance  split  levels  are  in  near  equilibrium  with 
each  other  up  to  about  85  km,  increasing  k^  (Mod3)  has  little  effect  on  the 
emission  in  the  region  where  good  data  are  found.  Above  85  km  the  reduction 
in  the  predicted  radiance  is  less  than  the  uncertainty  in  the  H20  contribution. 
For  these  reasons,  it  is  impossible  to  decide  from  this  model  whether  a 
larger  value  for  k^  is  appropriate. 

The  effect  of  reducing  k^  is  to  weaken  the  coupling  between  the  principal 
radiating  levels  and  the  side  levels,  and  thus  to  increase  the  populations 
of  the  former.  This  results  in  greater  emission  from  them.  By  reasoning 
similar  to  that  pertaining  to  k^*  we  conclude  that  the  value  for  k^  found 
in  the  basic  model  represents  a  lower  limit.  In  contrast,  larger  values  can¬ 
not  be  ruled  out.  An  order-of-magnitude  increase  in  k^  causes  a  decrease 
in  emission  which  removes  the  overall  H20-C02  prediction  from  agreement  with 

the  experiment.  However  a  larger  water  mixing  ratio - 2  ppm  in  the  50-70  km 

range,  which  is  still  consistent  with  the  SPIRE  results  in  the  6.3ym  region 

<15> - could  make  up  the  difference.  A  smaller  value  of  k^  might  have  the 

same  effect. 

One  result  which  has  not  been  mentioned  can  be  seen  in  Figures  7  and 
8:  in  a  range  of  altitudes  between  42  and  52  km  the  net  effect  of  collisions 

between  N2  and  C02  molecules  is  to  excite,  rather  than  deplete,  the  lOOlx 
levels  of  C02.  Generally  this  mechanism  quenches  high-lying  C02  levels.  In 
this  altitude  range  LTE  is  believed  to  prevail  for  the  low-lying  levels  in¬ 
volved  in  these  transitions  (e.g.,  lOOOx,  02201),  and  the  resulting  popula¬ 
tions  happen  to  be  great  enough  that  reactions  4  and  6  in  Table  2  run  in 


reverse. 


One  implication  of  this  result  is  that  for  paths  intersecting  this 
altitude  range  some  2. 7pm  emission  should  appear  even  for  night-time  con¬ 
ditions.  The  SPIRE  detectors  did  not  observe  night-time  emission  at  this 
wavelength.  This  is  explained  by  the  fact  that  in  the  absence  of  solar 
pumping  our  model  predicts  limb  radiances  from  these  bands  on  the  order  of 
10~9  watt/cm2-ster ,  which  is  well  below  the  noise  level  of  the  instruments. 

9 .  CONCLUSION 

A  model  for  the  prediction  of  C02  emission  in  the  2.7um  region  has  been 
developed,  and  tested  by  simulating  the  conditions  of  the  SPIRE  experiment. 

A  new  line-by-line  code  has  been  used  to  evaluate  the  rates  of  absorption 
of  solar  flux  in  many  relevant  bands.  These  rates  have  been  combined  with 
assumed  rates  for  inter-  and  intra-molecular  collisional  interactions  to 
evaluate  the  populations  of  high-lying  levels  responsible  for  emission  in 
the  101  combination  bands.  The  limb  radiance  has  been  calculated  from  these 
populations.  In  combination  with  emission  from  water  vapor  predicted  by 
another  model  <15>,  it  has  been  compared  with  the  SPIRE  dataset  for  daylight 
conditions.  Rate  constants  for  the  intramolecular  collisional  interactions 
have  also  been  altered  to  test  the  effect  of  possible  errors  in  their  values. 

In  general  the  model  has  been  found  to  give  good  agreement  with  the 
experiment.  It  has  also  established  the  need  to  include  the  intramolecular 
collisional  interactions  in  calculating  the  vibrational  populations  of  the 
high-lying  levels.  Among  our  specific  conclusions  are  that  (1)  the  resonant 
transitions  1001x-00001  dominate  the  limb  radiance  above  75  km;  (2)  the  flu¬ 
orescent  transitions  2001v-1000x  are  most  significant  on  lower  viewing  paths, 
but  several  other  sets  of  bands,  including  those  of  minor  isotopes,  contribute 
a  sizable  portion  of  the  observed  radiance;  (3)  the  principal  excitation  mecha¬ 
nism  is  solar  pumping  over  the  entire  range  of  altitudes  considered,  even 
when  the  line  centers  are  so  completely  bleached  out  that  absorption  occurs 
primarily  in  the  Voigt  wings;  (4)  the  high-lying  levels  are  de-excited 
primarily  through  radiative  processes  above  75  km  and  through  collisions 


with  N2  below  75  km;  (5)  the  intramolecular  collisional  transfer  rates  strongly 
affect  the  (integrated)  limb  radiance  and  also  the  spectral  shape  of  the  emis¬ 
sion;  and  (6)  the  values  of  the  rate  constants  used  to  characterize  these  in¬ 
teractions  could  be  adjusted  upward  without  causing  conflict  with  the  data, 
but  a  significant  downward  revision  appears  to  be  inappropriate. 
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III.  IMPROVEMENTS  IN  THE  NLTE  ALGORITHM 


1.  INTRODUCTION 

The  radiance  calculation  involved  in  the  comprehensive  CO^  model  and  in 
other  projects  which  we  have  undertaken  is  carried  out  by  program  NLTE.  In 
this  section  we  describe  modifications  which  were  recently  made  to  improve  the 
efficiency  and  running  speed  of  the  code. 

To  use  NLTE,  it  is  necessary  to  specify  a  band,  a  branch,  or  a  single  ro- 
vibrational  transition  of  a  particular  molecule  and  isotope.  One  then  provides 
an  atmospheric  profile  including,  most  importantly,  the  kinetic  temperature  and 
the  number  densities  of  the  vibrational  states  involved  in  the  transitions  (or 
equivalent  information) ,  and  also  data  from  the  A.F.G.L.  Atmospheric  Absorption 
Line  Parameters  Compilation  <l-3>  which  pertains  to  the  transitions  in  question. 
The  program  then  calculates  the  band  radiance  and  the  spectral  radiance  which 
a  hypothetical  sensor  would  detect  in  either  limb-  or  zenith-looking  geometry. 
Radiance  along  many  viewing  paths  may  be  determined  simultaneously. 

The  algorithm  upon  which  the  code  is  based  was  reported  in  1983  <4>,  along 
with  the  initial  release  of  the  program  and  instructions  for  its  use.  Augmented 
versions  of  this  first  release  have  been  applied  to  many  radiance  problems  in 
the  intervening  period.  The  relatively  long  execution  time  required  by  this 
detailed  calculation  has  been  only  a  minor  drawback  to  its  use.  In  the  future, 
however,  we  contemplate  its  use  in  multiple  or  iterative  applications,  within 
larger  programs.  This  makes  it  essential  to  have  more  efficient  code.  The 
changes  we  report  in  this  section  have  cut  the  run  time,  without  a  sacrifice  in 
accuracy,  by  factors  ranging  between  5  and  30  for  different  problems  we  have 
considered.  CO^  bands  can  now  be  run  in  less  than  20  seconds,  using  80  atmos¬ 
pheric  layers  and  5  viewing  paths. 


2.  PROGRAM  ORGANIZATION 


The  flow  chart  in  Figure  24  blocks  out  the  subtasks  performed  by  the  main 
program  after  the  input  and  initialization  steps  are  complete.  This  chart,  and 
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the  description  of  the  radiance  calculation  given  below,  are  meant  to  be  an 
introduction  only  to  the  organization  of  the  program  and  are  used  to  provide 
a  context  for  our  description  of  the  program  modifications.  The  physical  as¬ 
sumptions  and  the  quantities  we  actually  calculate  are  given  in  more  detail 
in  Reference  4. 

The  hierarchy  of  information  required  to  complete  the  radiance  calculation 
can  be  identified  as  follows:  (1)  quantities  which  are  independent  of  the 
transition  in  question,  such  as  atmospheric  profile  arrays;  (2)  line-dependent 
quantities  which  also  depend  on  altitude,  such  as  the  Doppler  width  and  corrected 
line  strengths;  (3)  line-  and  altitude-dependent  quantities  which  also  depend 
on  the  path,  such  as  path  column  densities;  and,  (4)  quantities  which,  in 
addition  to  the  above,  also  depend  on  the  frequency,  v.  Given  that  all  this 
information  is  available,  the  program  sets  V  equal  to  the  resonance  frequency, 
v  ,  and  begins  by  determining  the  radiance  at  the  near  edge  of  the  farthest 
layer  of  the  atmosphere  on  all  the  line-of-sight  paths. (See  Figure  25.)  Pro¬ 
ceeding  inward  toward  the  observer,  the  radiance  on  the  near  edge  of  the  second 
layer  is  the  sum  of  the  (attenuated)  contribution  from  the  first  layer  and  the 
contribution  from  the  second  layer  for  each  path.  This  summation  continues 
from  layer  to  layer  until  the  contribution  from  the  last  (nearest)  layer  is 

added  in  to  give  the  observed  radiance  (in  photons/sec-cm2-ster-cm-1)  at  V  . 

o 

These  steps  constitute  block  D  in  Figure  24.  Then  v  is  changed  and  the  procedure 
repeats.  When  a  sufficient  number  of  points  have  been  evaluated  to  provide  a 
complete  radiance  profile,  or  spectral  lineshape,  a  numerical  integration  over 
frequency  is  performed  to  get  the  total  line  radiance.  In  the  end,  the  con¬ 
tributions  from  all  the  lines  give  the  band  radiance. 

One  can  see  that  the  radiance  calculation,  block  D,  with  its  loops  over 
altitude  and  path,  contains  the  most  deeply  embedded  calculations  the  program. 

It  turns  out,  however,  that  the  most  time-consuming  calculations  are  those  of 
the  Voigt  lineshape  function,  block  C,  which  is  also  deeply  embedded  but  which 
loops  only  over  altitude.  The  obvious  strategy  for  improving  the  speed  of  the 
program  is  therefore  to  substitute  a  faster  Voigt  profile  routine  and  to  reduce 
the  number  of  calls  to  this  routine.  These  tactics  in  fact  are  responsible  for 


the  enhanced  program  efficiency  we  are  reporting.  In  Section  3  we  describe  our 
adaptation  of  an  existing  Voigt  profile  algorithm  which  is  five  to  six  times 
faster,  on  average,  than  the  procedure  we  used  previously.  In  Section  A  we 
describe  a  new  semi-adaptive  integration  algorithm  which  requires  far  fewer 
points  than  the  compounded  Simpson's  rule  previously  in  use.  We  also  describe 
an  analytical  approximation  to  the  tail  of  the  radiance  profile  which  eliminates 
the  need  to  carry  the  computation  very  far  from  the  line  center. 

3.  VOIGT  PROFILE  ALGORITHMS 
3.1.  Lineshape  Considerations 

One  physical  assumption  of  the  NLTE  algorithm  is  that  the  frequency-dependence 

of  the  absorption  cross-section  (and  hence  the  emission  lineshape)  is  given  by  a 

Voigt  profile,  which  is  the  convolution  of  the  Doppler  lineshape  (a  Gaussian) 

and  the  pressure-broadened  lineshape  (a  Lorentz  function).  At  altitudes  where 

non-LTE  conditions  prevail  the  Lorentz  component  is  often  very  weak,  because  of 

the  relative  infrequency  of  collisions  in  the  rarified  atmosphere.  Even  with  a 

-2  -3 

Lorentz-to-Doppler  linewidth  ratio  of  10  or  10  ,  however,  we  have  shown  that 

use  of  the  Voigt  profile,  rather  than  the  Doppler  lineshape  which  is  easier  to 
compute,  is  extremely  important  because  of  the  much  slower  dropoff  in  the  wing 
of  the  Voigt  function.  That  is,  many  lines  may  be  severely  self-absorbed  in  the 
line-center  even  for  viewing  paths  traversing  only  very  high-altitude  regions; 
in  such  cases  the  contribution  in  the  wings,  where  the  self-absorption  is  much 
less,  constitutes  a  large  part  of  the  total  radiance  profile,  and  it  is  necessary 
to  have  an  accurate  representation  of  the  lineshape  in  the  wings.  Errors  of  501 
to  100%  can  easily  be  made  by  neglecting  the  collision-broadening. 


3.2  Basic  Representation 


The  Lorentz,  Doppler,  and  Voigt  profiles  are  defined  by 
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V(v)  =  L(v)  ©  D(v)  (2) 

where  a  and  a  are  the  pressure-broadened  and  Doppler  linewidths,  re- 
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spectively,  and  is  the  resonance  frequency.  The  linewidths  are  altitude- 
dependent  through  the  pressure  and  temperature  <1>.  Since  L  and  D  are  each 
normalized  to  unity,  V  is  also  normalized.  By  making  the  substitutions 
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the  following  form  is  achieved: 
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V(v)  is  thus  parameterized  by  dimensionless  quantities,  y  and  x,  representing 
the  linewidth  ratio  and  the  frequency.  It  happens  to  be  proportional  to  the 
real  part  of  a  complex  function,  w(z). 
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=  —  / -  Re[w(x  +  iy)] 
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where  w(z)  is  related  to  the  complex  error  function  by  <5> 

e  2  (l  -  erf  (-iz)j  (6) 

and  x  and  y  are  given  above.  Many  algorithms  for  the  Voigt  function  are  based 
on  this  representation:  one  finds  an  efficient  way  of  calculating  Re[w(z)] 
starting  with  Eq.  (6). 


— OO 


3. 3  Rybicki's  Algorithm 

Rybicki's  algorithm  <6>,  which  was  used  in  the  first  release  of  NLTE  as 
funct  ~i  VWERF,  determines  Re[w(z)]  by  expanding  e  in  a  sampling  representation 
and  evaluating  the  resulting  integrals  by  the  method  of  residues.  This  gives  a 
series  which  is  cut  off  at  31  terms.  It  is  a  highly  accurate  representation, 
giving  between  8  and  12  correct  significant  figures  over  the  full  range  of 
parameters  we  tested.  Since  our  accuracy  requirements  are  not  nearly  so  stringent 
as  this,  we  tried  cutting  off  the  series  with  fewer  terms.  The  speed  is  worse 
than  directly  proportional  to  the  number  of  terms,  however,  and  the  great  accuracy 
quickly  disappears  for  certain  portions  of  the  x-y  plane  when  fewer  terms  are  used. 
A  19-term  series  is  still  advantageous,  but  other  methods  proved  to  be  better. 

Rybicki's  algorithm  makes  use  of  an  asymptotic  formula  for  very  large  x. 


3.4  Pierluissi's  Algorithm 

Pierluissi  et  al  <7>  developed  a  fast  Voigt  algorithm  which  we  have  modified 
and  adapted  to  our  purposes. 
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Pierluissi's  algorithm,  realized  as  function  ZV01GT,  divides  the  x-y  plane 
into  three  regions.  Region  I  is  specified  by  (0<x<3,  0<y<1.8).  Region  il  is 
bounded  by  region  I  on  one  side  and  by  x=5,  y=5  on  the  other.  Region  111  is 
specified  by  (x>5,  yN5)  . 

In  Region  I,  Re[w(z)]  is  determined  by  using  a  series  expansion  <5>  for  the 

complex  error  function  in  Eq.  (6) .  The  number  of  terms  in  the  series  is  a 

linear  function  of  x.  The  maximum  number  of  terms  is  29.  Each  call  requires 

evaluations  of  one  (real)  sine,  cosine,  and  exponential.  The  series  computation 

is  between  two  and  three  times  faster  than  VWERF  in  the  comparable  x-y  region. 

Its  accuracy  is  practically  independent  of  y  for  y<l ,  but  is  very  sensitive  to 

the  value  of  x.  In  fact,  for  small  x,  extremely  accurate  values  are  obtained. 

-3  -3 

Fractional  errors  of  10  and  5x10  are  obtained  for  x<3.6  and  x<3.8,  respec¬ 
tively.  It  is  not  useful  for  y£2. 

In  Regions  II  and  III  the  evaluation  of  Rc[w(z)]  is  accomplished  by  6-  and 
4-point  Gauss-Hermite  quadrature  <8,9>.  Since  the  weighting  function  for  this 
quadrature  is  exactly  the  exponential  appearing  in  the  integral  in  Eq.  (6),  no 
transcendental  function  evaluations  are  required  at  all  and  the  calculation  is 
very  fast. 

Generally  speaking,  the  accuracy  of  the  4-  and  6-point  rules  is  excellent 
for  x>5.  It  improves  with  increasing  x  and,  for  x>6,  does  not  depend  on  y  for 
tne  range  of  y  values  tested.  It  rapidly  worsens  with  decreasing  values  of  x, 
however,  especially  for  small  values  of  y.  In  fact,  within  Region  II  the  6-point 
rule  gives  fractional  errors  of  up  to  10%  for  y=.006  and  practically  100%  (e.g., 
completely  wrong  results)  for  values  of  .001  and  less.  It  is  conceivable  that 
this  fact  was  overlooked  in  the  original  formulation  because  of  the  presumed 
unimportance  of  the  wings  for  small  y-values.  However,  many  strong  infrared 
bands  contain  lines  for  which  y  is  less  than  .001  at  altitudes  for  which  self¬ 
absorption  is  very  important.  The  need  to  accurately  compute  the  wing  contribu¬ 
tions  therefore  made  it  necessary  for  us  to  modify  the  Pierluissi  algorithm  be¬ 
fore  using  it. 
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3.5  Tests  for  Accuracy  and  Speed 

In  order  to  test  different  Voi gt-function  algorithms,  we  wrote  a  program 
to  do  a  brute-force  integration  of  Eq.  (4)  to  serve  as  a  benchmark,  and  chose 
a  range  of  parameters  within  which  to  evaluate  the  functions.  The  ranges  are 
. 000001 ^yS10  and  05x<30. 

All  tests  were  performed  on  a  CYBER  750  machine  running  under  the  NOS 
operating  system. 

Figure  26  gives,  as  a  function  of  altitude,  an  approximate  range  of  values 
of  the  parameter  y  which  are  encountered  when  emission  in  two  of  the  strongest 
CO^  bands  is  calculated.  That  is,  for  example,  for  a  given  altitude  any  line 
in  the  various  bands  will  be  characterized  by  a  Lorentz  to  Doppler  linewidth 
ratio  falling  between  the  two  top  curves  in  Figure  26.  From  this  and  the  view¬ 
ing  path  one  can  readily  determine  the  range  of  y-values  for  which  the  Voigt 
profile  will  need  to  be  evaluated  for  any  particular  run.  We  use  the  U.S. 
Standard  Atmosphere  1976  <10>  for  our  pressure-versus-al titude  profile. 

We  made  some  test  runs  with  NLTE  with  the  CO^  v,  fundamental  band  at  15lim 
in  limb  geometry  and  counted  the  number  of  calls  to  the  Voigt  function  act  rding 
to  ranges  of  values  of  x  and  y.  The  purpose  of  this  exercise  was  to  make  it 
possible  to  optimize  our  algorithm  for  speed  in  the  regions  of  the  x-y  plane 
encountered  most  frequently  in  typical  runs.  Table  4  gives  the  percentages  of 
calls  in  different  regions  for  a  run  in  which  paths  with  tangent  heights  of  70, 
80,  90,  and  100  km  were  done  simultaneously.  Table  5  gives  the  percentage  of 
calls  within  different  ranges  of  the  parameter  x  alone,  for  jobs  in  which  single 
paths  with  tangent  heights  between  50  and  100  km  were  run.  In  all  cases,  the 
atmosphere  was  divided  into  one-kilometer  layers  between  the  tangent  point  and 
an  altitude  of  150  km.  For  each  viewing  path,  lines  for  which  the  analytic 
"thin-line"  approximation  is  correct  to  1%  or  better  were  disregarded  in  these 
tests . 

The  results  in  Tables  4  and  5  are  predictable  in  several  respects,  notably 
that  the  lower  viewing  paths  reauire  a  great  many  calls  in  the  region  x'5.  This 
occurs  because  of  the  very  slow  dropoff  of  the  Voigt  function  in  the  far  wings: 
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Altitude  (km) 


Minimum  and  maximum  values  of  the  Voigt  parameter, 
y,  as  a  function  of  altitude,  for  the  lines  in  the 
CO?  v 2  bands  near  1 5um  and  for  the  lines  in  the  CO 
v i +v 3  bands  near  2.7pm.  Values  for  the  CO?  v?  bam 
not  shown,  are  slightly  greater  than  those  of  the 
Vi+O-j  bands. 


\  X 

Percentage 

y  \ 

0-1.67 

1.67-3.33  3. 

33-5.00 

>5 

io~3  -  10"2 

4.8 

3.6 

2.6 

4.9 

10'4  -  10“ 3 

4.8 

3.7 

2.6 

5.4 

10" 5  -  io~4 

4.9 

3.7 

2.6 

5.0 

io-6  -  10-5 

7.3 

5.2 

3.1 

4.7 

0  -  10*6 

14.2 

8.2 

4.1 

4.6 

Total 

36.0 

24.4 

15.0 

24.6 

Table  #4  Percentage  of  calls  to  the  Voigt  function  for  different  ranges  of 

x  and  y,  for  a  run  using  the  V2  fundamental  of  CO2  in  limb  geometry. 
Paths  with  tangent  heights  of  70,  80,  90,  and  100  km.  were  evaluated 
simultaneously.  The  total  number  of  calls  was  286,000. 


Percentage 


hT\ 

0-1.67 

1.67-3.33 

3.33-5.00 

>5 

N 

(km) 

50 

16.8 

13.0 

10.7 

59.5 

889,700 

60 

20.9 

15.2 

11.5 

52.4 

596,970 

70 

43.1 

28.2 

14.8 

13.9 

239,760 

80 

61.6 

31.3 

7.1 

0 

136,150 

90 

67.5 

27.9 

4.6 

0 

102,540 

100 

72.7 

24.7 

2.6 

0 

77,250 

Table  // 5  Percentage  of  calls,  and  total  number  (N)  of  calls,  to  the  Voigt 

function  for  four  ranges  of  x,  for  runs  using  the  v2  fundamental 
of  CO2  in  limb  geometry  and  different  tangent  heights,  h  . 
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for  the  very  thick  lines  the  optical  depth  as  far  out  as  25  Doppler  widths  may 

-2 

He  greater  than  unity  even  with  y-values  on  the  order  of  10  .  Conversely, 

for  higher  paths  the  contributions  in  the  far  wings  are  much  less  significant 
ana  the  program  can  cut  off  the  numerical  integration  much  closer  to  the  line 
center.  These  results  suggest  that  the  design  of  an  algorithm  for  use  in  a 
multitude  of  such  applications  should  allow  for  fast  execution  both  near  the 
line  center  and  in  the  far  wings,  perhaps  at  the  expense  of  speed  in  the  region 
2<x<5. 

We  tested  many  Voigt  algorithms  for  accuracy.  These  included  the  three 

separate  procedures  embodied  in  ZVOIGT  <7>,  a  2-point  Gauss-Hermite  quadrature 

formula  < 1 1 > ,  and  the  substitution  of  a  simple  Doppler  profile.  None  of  these 

provides  an  acceptable  representation  of  the  Voigt  profile  over  the  entire 

plane,  or  even  over  as  much  of  it  as  Rybicki's  sampling  representation  does  in 

VWERF.  In  the  end  we  pieced  together  several  procedures  satisfying  the  criterion 

-3 

that  the  Voigt  profile  be  calculated  with  a  fractional  error  less  than  10  for 

-3 

0Sx'il.2  and  less  than  5x10  for  x>1.2. 

The  Gauss-Hermite  formulas  are  very  fast.  They  are  accurate  for  large  :: 

(more  or  less  independent  of  y)  and  also  for  large  y.  For  a  given  parameter 

set  (x,y)  the  6-point  rule  is  of  course  more  accurate  than  the  4-  or  2-point 

rules.  Figures  27-29  show  the  regions  of  the  x-y  plane  where  these  procedures 

produce  values  more  accurate  than  a  certain  threshold  value,  F.  We  have  con- 

-4  -3  -3 

sidered  fractional  errors  of  10  ,10  ,  and  5x10  for  F.  (Actually,  the  true 

boundaries  are  not  nearly  so  smooth  as  we  have  drawn  them  in  these  figures. 
Because  the  quadrature  formulas  are  polynomial  fits  to  the  Voigt  function,  there 
are  pockets  and  islands  in  the  x-y  plane  to  the  left  of  the  boundaries  within 
which  the  function  is  also  represented  accurately.  The  accuracy  criteria  are 
satisfied,  however,  everywhere  to  the  right  of,  and  above,  these  curves.) 

Not  surprisingly,  the  Doppler  profile  D(v)  is  a  useful  substitution  for 
V(v)  for  small  values  of  x  and  y.  Figure  30  gives  the  "regions  of  validity" 
for  this  substitution  for  the  same  three  threshold  values. 
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Figure  27.  Regions  of  the  x-y  plane,  within  which  a  2-point 
Gauss-Hermite  quadrature  routine  approximates  the 
Voigt  function  with  a  fractional  error  less  than 
F,  lie  above  and  to  the  right  of  the  three  curves 
shown.  The  curves  correspond  to  F  =  10—  4 ,  10~3, 
and  5xlO~3. 
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Figure  28. 


Regions  of  the  x-y  plane,  within  which  a  4-point 
Causs-Hermi te  quadrature  routine  approximates  the 
Voigt  function  with  a  fractional  error  less  than 
F,  lie  above  and  to  the  right  of  the  three  curves 
shown.  The  curves  correspond  to  F  =  10“3, 

and  5x10“ 3 . 
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Figure  29.  Regions  of  the  x-y  plane,  within  which  a  6-point 
Gauss-Hermite  quadrature  routine  approximates  the 
Voigt  function  with  a  fractional  error  less  than 
F,  lie  above  and  to  the  right  of  the  three  curves 
shown.  The  curves  correspond  to  F  =  10-4,  10-3, 
and  5x10“ 3 


Figure  30.  Regions  of  the  x-y  plane,  within  which  the  Doppler 
profile  can  be  substituted  for  the  Voigt  function 
with  a  fractional  error  less  than  F,  lie  below  and 
to  the  left  of  the  three  curves  shown.  The  curves 


correspond  to  F  =  10" 


10 


-3 


and  5x10 
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The  ZVOIGT  series  used  in  Pierluissi's  region  I  is,  as  noted  above,  very 
accurate  for  small  x.  Its  accuracy  deteriorates  rapidly  with  increasing  x, 
and  it  is  not  good  enough  for  the  series  to  be  used  at  all  for  y>2.  In  order 
to  extend  its  useful  range  into  the  region  (x>3,  y<001)  we  carried  it  out  to  a 
maximum  of  55  terms  and  chose  a  linear  function  of  x  to  determine  the  number 
of  terms  to  be  taken  which  is  somewhat  different  from  Pierluissi's.  With  its 
three  evaluations  of  transcendental  functions  and  a  multitude  of  arithmetic 
operations,  this  procedure  is  slower  than  the  others  mentioned  above. 

Table  6  gives  the  average  execution  time  for  six  FORTRAN  routines  used  to 
compute  V(v)  on  the  CYBER  750.  The  routines  were  compiled  in  the  fast-execute 
mode  (0PT=2  on  FTN5) .  GH2,  GHA,  and  GH6  are  the  Gauss-Hermite  routines.  DOPP 
evaluates  the  Gaussian  function,  D(v) .  SERIES  is  the  ZVOIGT  series,  modified 
as  described  above.  It  is  the  only  one  of  these  whose  execution  time  depends 
on  the  parameters,  as  the  series  has  fewer  terms  when  x  is  small.  VWERF  is 
Rybicki's  algorithm. 

In  each  case  cited  above  the  only  overhead  operations  are  the  single  call 
and  return.  That  is,  in  the  timing  tests  there  were  no  branching  operations  or 
extra  arithmetic.  Table  6  shows,  with  the  entries  DUMMY  and  GO  TO,  respectively, 
the  approximate  time  for  the  call  and  return,  and  the  approximate  extra  time 
used  by  a  conditional  GO  TO  statement.  Also,  a  real  transcendental  function 
takes  approximately  10  microseconds  to  execute  within  any  procedure. 

3.6  Function  ZVGT 

Our  Voigt  algorithm,  realized  as  function  ZVGT,  is  a  combination  of  four 
separate  procedures:  the  2-  and  6-point  Gauss-Hermite  formulas,  the  Doppler 
substitution,  and  the  extended  ZVOIGT  series.  Figure  31  indicates  the  particular 
division  of  the  x-y  plane  we  use  in  this  procedure.  The  division  is  implemented 
by  several  conditional  GO  TO  statements.  One  may  note  that  the  Doppler  and  GH6 
regions  could  be  expanded  into  portions  of  the  region  where  the  slower  series 
operates.  The  cost  of  doing  this,  however,  is  in  each  case  the  evaluation  of  a 
logarithmic  or  exponential  function  before  effecting  the  conditional  branches. 


Procedure 


Time 

(ysec) 


GH2  13 
GH4  15 
GH6  17 
DOPP  19 
SERIES  97 
VWERF  230 


DUMMY  9 


GO  TO  2.5 

TRAN SC  10 


Table  //6  Average  execution  times,  in  microseconds,  for  various  FORTRAN 
procedures  and  operations  on  the  CYBER  750.  See  text  for  a 
full  explanation. 


N.  x 

y  ^ 

Average 

0-1.67 

Execution 

1.67-3.33 

Time  (ysec),  ZVGT 

3.33-5.00 

>5 

1.5 

-  10 

21 

21 

21 

15 

.1 

-  1.5 

70 

98 

123 

16 

10~2 

-  io-1 

71 

94 

121 

15 

10~3 

CM 

1 

O 

1 

70 

97 

122 

15 

io-4 

-  10  3 

23 

92 

120 

16 

io“5 

-  io-4 

23 

94 

122 

15 

io-6 

1 

C 

*— < 

1 

23 

93 

122 

16 

0 

-  io'6 

23 

23 

23 

15 

Table  // 7  Average  execution  times  (ysec)  for  function  ZVGT  for  different 
regions  of  the  x-y  plane. 
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Figure  31.  Partition  of  the  x-y  plan  into  regions  for  function  ZVGT. 

2-point  Gauss-Hermite  quadrature  is  used  for  x>4.9.  The 
6-point  rule  is  used  for  (y>1.5,  x<4.9).  The  Doppler 
function  is  used  for  (x<1.7,  y<.0009)  and  for  (y<10-6. 


The  extra  time  required  for  these  functions  is  greater,  in  a  typical  run  with 
parameters  distributed  as  in  Table  A,  than  the  time  saved  by  implementing  the 
faster  arithmetic. 

As  noted  above,  ZVGT  calculates  V(v)  with  a  fractional  error  of  less  than 
-3  -3 

10  for  x$1.2  and  5x10  for  x>1.2.  That  is,  the  function  is  accurate  to  at 
least  .1%  in  the  line  center  and  .5%  elsewhere.  In  fact,  these  limits  are 
approached  only  at  the  boundaries  of  the  sections.  Over  most  of  the  x-y  plane 
the  routine  is  much  more  accurate  than  this,  especially  for  small  x. 

Table  7  gives  the  average  execution  time  required  by  a  single  call  and 
return  to  function  ZVGT  for  different  regions  of  the  x-y  plane.  For  the  distri¬ 
bution  of  calls  given  in  Table  A,  the  average  execution  time  is  about  A6  usee. 

For  the  distribution  of  calls  required  for  tangent  heights  between  50  and  100  km, 
the  average  execution  time  ranges  between  36  and  A9  usee  using  ZVGT.  Because  of 
the  emphasis  on  very  large  x,  where  the  routine  is  very  fast,  the  lower  tangent 
paths  require  a  smaller  average  time;  of  course,  they  also  require  a  larger 
number  of  calls. 

We  note  that  by  calculating  a  large  array  of  values  at  one  time  rather  than 
calling  the  routine  for  each  sepeate  evaluation,  an  additional  9  ysec  could  be 
saved,  per  average  evaluation. 

The  source  listing  of  ZVGT  is  given  in  Appendix  A. 

A .  NUMERICAL  INTEGRATION  ALGORITHM 
A  .  1  General  Description 

2 

Nl.TE  computes  the  band  radiance,  in  watt/cm  -ster,  bv  performing  a  numerical 
integration  over  the  spectral  radiance  profile  computed  for  each  line,  and  then 
summing  over  the  lines  in  the  band.  In  the  first  release,  the  numerical  inte¬ 
gration  procedure  was  awkward  to  use  because  it  was  necessary  to  specify  both  a 
frequency  interval  sufficiently  small  for  the  compounded  Simpson's  rule  to  be 
accurate,  and  also  a  number  of  integration  points  sufficiently  large  to  avoid  a 
significant  truncation  error  in  the  far  wing  of  the  strongest  lines.  This  in¬ 
evitably  led  to  unnecessary  work  in  the  wings  of  weak  lines. 


I 

I 

I 

We  have  replaced  the  original  integration  algorithm  with  a  new  procedure  j 

which  avoids  these  input  requirements  and  also  is  more  efficient  for  three 

I 

reasons.  First,  it  uses  Gauss-Legendre  integration  formulas,  which  are  in¬ 
herently  more  accurate  than  Simpson's  rule,  given  a  certain  number  of  points  I 

% 

in  a  fixed  interval.  This  allows  us  to  use  fewer  points  throughout  the  J 

profile.  Second,  it  is  sensitive  to  the  characteristics  of  individual  lines 

rather  than  to  those  of  an  entire  band,  so  it  cuts  off  the  calculation  at 

distances  from  the  line  center  which  are  correct  for  both  strong  lines  and 

weak  lines.  Moreover,  it  recognizes  that  numerical  integration  in  the  far 

wings  of  the  radiance  profile,  where  the  integrand  decreases  monotonically 

and  quite  slowly,  can  be  carried  out  with  fewer  points  per  frequency  interval 

than  it  can  nearer  the  center  where  the  function  may  have  a  non-central  maximum 

and  in  any  event  varies  more  rapidly.  This  allows  us  to  use  still  fewer  points 

per  frequency  interval  in  the  wings.  Third,  we  use  an  analytical  approximation 

to  evaluate  the  integrated  contribution  of  the  far  tail  of  the  radiance  profile, 

thereby  eliminating  the  truncation  error  and  allowing  the  numerical  evaluation 

to  be  cut  off  closer  to  the  line  center. 

On  the  whole,  these  three  changes  have  succeeded  in  inducing  the  number  of 
function  evaluations,  and  hence  the  run  time,  by  a  considerable  amount.  They 
are  discussed,  one  by  one,  in  the  following  sections.  Some  of  the  work  described 
below  relies  heavily  on  analysis  first  presented  by  Bullitt  et  al  <12>. 

4 . 2  Integration  Rule 

We  use  Gauss-Legendre  formulas  for  our  numerical  integration  in  NLTE  because 
of  their  high  accuracy  <8>.  Two-,  four-,  or  eight-point  rules  can  be  specified. 
Built-in  flexibility  is  provided  because  the  program  divides  the  frequency  range 
into  successive  "panels",  or  intervals,  within  which  it  applies  the  Gauss  formula. 

The  first  panel  comprises  the  interval  from  the  line  center  out  to  one  "standard 
Doppler  width".  For  all  paths,  the  integrated  radiance  in  this  interval  is 
determined  before  the  second  panel  is  considered.  The  second  panel  comprises  the 
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Fig.  32.  Radiance  profiles  of  4  lines  of  the  R  branch  of  the 

C02(v7)  fundamental  near  667  cm-1.  The  curves,  which 
were  calculated  by  program  NLTE  were  plotted  for  limb- 
look  geometry  f or  a  tangent-height  of  70kra.  The  solid 
lines  were  found  assuming  a  Voigt  lineshape.  The  dotted 
line  is  the  R2  profile  assuming  a  Doppler  lineshape. 

The  arrow  at  the  bottom  indicates  the  Doppler  linewidth 
at  the  tangent-height. 
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frequency  interval  between  one  and  two  "standard  Doppler  widths",  and  so  on. 

This  gives  us  an  accumulating  total  for  the  integrated  intensity  of  the  line 
as  we  proceed,  which  is  useful  for  evaluating  the  cutoff  criterion  in  the  wing. 

A  "standard  Doppler  width"  is  defined  as  a  Doppler  width  determined  by 
using  the  band  center  for  the  frequency  and  the  lowest  temperature  in  the  alti¬ 
tude  range  under  consideration  for  the  kinetic  temperature. 

We  have  found  that  the  default  value  of  four  points  per  panel  is  sufficient 
for  very  accurate  results.  Use  of  two  points  per  Doppler  width  often  gives 
errors  in  the  band  radiance  which  are  less  than  errors  introduced  by  other 
approximations  built  into  the  program. 

A . 3  Adjustments  for  Different  Lines 

Figure  32  illustrates  the  great  differences  which  characterize  radiance 
profiles  of  different  lines  within  the  same  band,  provided  some  of  the  lines 
are  self-absorbed  along  the  viewing  path.  The  weakest  lines  mimic  Doppler 
lines;  stronger  lines  decrease  monotonically  from  the  center  but  have  much  more 
pronounced  wings;  and,  for  typical  temperature  profiles,  the  strongest  lines 
have  maxima  and  possibly  minima  away  from  the  line  center,  as  well  as  wings 
where  the  function  decreases  smoothly  and  monotonically  but  very  slowly  over  a 
large  range  of  frequencies. 

As  mentioned  above,  four  points  per  Doppler  width  usually  give  accurate 
results  for  the  numerical  integration  over  the  whole  profile.  Since  the  regions 
where  the  integrand  varies  most  rapidly,  which  in  general  are  the  central  region 
and  wherever  else  extrema  occur,  are  the  most  difficult  to  integrate  accurately, 
it  stands  to  reason  that  the  regions  where  the  integrand  is  smoothest  require 
fewer  points  than  this.  In  our  case,  the  wings  of  the  radiance  profile  can  be 
integrated  very  efficiently  by  expanding  the  panel  width  and  using  the  same 
number  of  points  per  panel  as  near  the  line  center.  Of  course,  it  is  necessary 
to  be  able  to  specify  the  frequency,  beyond  which  the  radiance  profile  is  smooth 
enough  to  expand  the  panels,  on  a  line-by-line  basis. 


Bullitt  et  al  <12>  have  shown  that  the  radiance  profile  can  be  represented 
as  the  product  of  two  quantities.  The  first  is  the  absorptivity,  (1  -  expCi^)), 
where  T  is  the  optical  depth  along  the  line-of-sight  path  at  frequency  v.  The 
second  is  <R>/<Roo>,  where  <R>  is  an  effective  blackbodv  radiance  averaged  over 
the  path  with  a  certain  weighting  function  <12>,  and  <Ro0>  is  its  limiting  value 
in  the  far  wing. 

The  second  of  these  two  factors  smoothly  approaches  unity  in  the  wings,  so 
the  functional  dependence  of  the  radiance  profile  there  is  essentially  that  of 
the  absorptivity.  In  fact,  for  t  «1,  one  sees  by  expanding  the  exponential 
that  the  frequency-dependence  is  entirely  contained  in  T  ,  which  is  proportional 
to  V (v) . 

On  the  other  hand,  if  t  >>1  in  the  line  center,  the  absorptivity  is  nearly 
constant  (and  nearly  unity)  so  the  functional  dependence  of  the  radiance  profile 
is  basically  that  of  the  second  factor.  This  second  factor,  being  akin  to  a 
blackbodv  radiance,  is  influenced  by  the  temperature  along  the  path  and  is 
thereby  responsible  for  maxima  and  minima  appearing  near  the  center  of  the 
profile  if  ",'-•>1.  That  is,  the  observer  "sees"  emission  from  portions  of  the 
path  only  a  certain  distance  into  the  atmosphere  because  of  self-absorption, 
and  the  amount  of  radiation  reaching  him  depends  on  the  temperature  along  those 
portions.  Since  the  portions  "seen"  depend  on  frequency  through  the  optical 
depth,  the  radiance  profile  may  have  extrema  if  the  temperature  profile  varies 
rap i d 1 v . 

The  considerations  introduced  in  the  last  two  paragraphs,  which  are  illus¬ 
trated  quite  effectively  in  Figure  9  of  Reference  12,  imply  that  beyond  the 
frequency  at  which  the  absorptivity  departs  significantly  from  its  value  at  V  , 
there  will  be  no  extrema  in  the  radiance  profile.  That  is,  for  all  frequencies 
greater  than  some  .  such  that  T  1,  the  observer  will  "see"  the  entire  path  and 
the  f  requency-dependence  of  -R  /*-  R  will  he  slight;  however,  this  is  exactly 
the  range  in  which  the  monotonic  decrease  in  the  absorptivity  begins  to  dominate 
the  function.  For  these  reasons  the  program  uses  t  - . 5  as  a  flag  for  the 
frequency  beyond  which  the  functional  dependence  of  the  radiance  profile  is 
sufficiently  slowly-varying  t  ci  expand  the  panel  width.  This  is  actually  a 
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fairly  conservative  benchmark;  we  find  that  the  outermost  peak  in  typical 

thick-line  radiance  profiles  comes  between  i  - 2  and  t  -5. 

v  v 

The  rule  for  expanding  the  panel  width  was  established  by  trial  and  error. 

We  double  the  panel  width  for  each  panel  occurring  past  the  aforementioned 

"flag".  This  obviously  carries  the  numerical  integration  into  the  far  wings 

very  rapidly.  The  reason  it  works  as  well  as  it  does  can  be  seen  in  the  error 

term  for  four-point  Gauss-Legendre  quadrature,  which  for  a  fixed  interval  size 

goes  as  the  eighth  derivative  of  the  integrand  at  some  point  in  the  integration 

interval.  For  the  wings,  as  we  pointed  out  earlier,  the  integrand  goes  as 

V(v),  which  is  very  nearly  proportional  to  the  Lorentzian  L(v)  there.  The 

eighth  derivative  of  L(v)  is  proportional  to  (v-v  )  With  such  a  rapidly 

o 

decreasing  error  for  fixed  intervals,  one  can  see  how  the  intervals  can  be 
expanded  rapidly  as  ( v-\ ;  )  increases  without  compromising  the  accuracy  of  the 
calculation. 

Weak  lines  with  T^<<1  at  the  line  center  can  be  evaluated  by  the  "thin-line" 
approximation  which  neglects  absorption  along  the  viewing  path.  When  Tv<.5  and 
the  thin-line  approximation  is  not  used  the  panel  widths  are  not  doubled  until 
the  fourth  panel  is  reached.  Usually  the  truncation  procedure  is  invoked  quite 
quickly  for  such  lines,  anyway. 

4 . 4  Truncation  Procedure 

The  third  component  of  our  revised  procedure  for  integrating  the  radiance 
profile  is  a  cutoff  of  the  numerical  evaluation.  The  radiance  profile  is  of 
course  originally  written  as  an  integral  over  the  viewing  path  <12'.  If  s" 
parraeterizes  the  path  and  s  and  s'  are  its  near  and  far  ends,  respectively,  the 
radiance  is  <4,  12> 

A  s 

Iv(s)  =  ^  j  nu(s")Vv(s")  exp  (-T v (s ,s") ) ds "  (7) 

s' 

where  A  is  the  F.instein  coefficient  for  spontaneous  emission,  n  is  the  number 

3  u 

density  of  the  upper  radiating,  .state  in  molecules/cm  ,  and  the  path-dependence 


of  the  Voigt  profile  has  been  made  explicit.  At  the  path's  end,  the  integrated 
radiance  in  the  tail,  that  is  in  the  region  |v-\)q  |<°°  is 


T(vc,s)  =  2  j°°  Iv(s)dv 

v 

c 


(8) 


which,  with  T  «1 ,  can  be  written  as 

v 


T(v  ,s) 
c 


A_ 

2tt 


s 

n  (s") 
u 

* 

s 


oo 

j  Vv(s")  dv 

V 

c 


ds” 


(9) 


Using  any  of  several  representations  of  the  tail^of  the  Voigt  function,  one  can 
show  that  the  v-integral  is  equal  to  — v~)  to  or<*er  (y/*c)^  where  y  and 

x  are  the  Voigt  parameters  of  Eqs.  (3)  and  x^  is  x  evaluated  at  v=vc*  This 
quantity  is  of  course  altitude-dependent  through  the  Lorentz  width,  a^,  and  must 
be  integrated  over  s".  In  NLTE,  the  s"  integral  is  just  a  sum  over  homogeneous 
layers,  however,  and  is  done  on  the  initial  pass  through  the  atmosphere  as  the 
thin-line  approximation  is  evaluated  for  each  path.  Given  this,  it  is  easy  to 
calculate  T(\'c,s)  for  any  given  for  all  the  paths. 

The  difficult  part  is  determining  how  accurate  the  truncation  term  is.  Our 
procedure  for  determining  where  to  cut  off  the  numerical  integration  and  invoke 
the  truncation  approximation  is  based  on  empirical  tests  and  is  likely  to  be 
refined  in  the  future.  At  present  we  evaluate  the  truncation  term  at  the  end  of 
each  panel  past  the  third  one  and  compare  it  with  the  quantity  a,  defined  as 


6ei 

c 

a  =  - 

i 

v 

c 

where  I  is  the  cumulative  integrated  radiance  in  the  panels  treated  so  far,  T 

c  v 

Q 

is  the  optical  depth  for  the  path  at  the  end  of  the  last  panel,  and  C  is  the 
desired  accuracy  of  the  integrated  radiance  in  the  line,  expressed  as  a  fractional 
error.  If  a  is  greater  than  Tfv^.s),  the  expression  for  the  tail  is  presumed  to 
be  valid  and  the  numerical  evaluation  of  the  radiance  profile  is  cut  off. 
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IV.  ATMOSPHERIC  AURORAL  RADIANCE  CODE 


1.  INTRODUCTION 

Over  the  past  several  months,  much  effort  has  been  devoted  to  the  construc¬ 
tion  of  the  computer  program,  AARC,  whose  overall  purpose  is  the  calculation  of 
infrared  emissions  from  nitric  oxide  excited  by  energetic  auroral  electrons  and 
background  earthshine  radiation. 

Many  of  the  procedures  involved  in  the  calculation  are  new,  while  others  had 
already  existed.  Several  improvements  and  corrections  were  made  to  these  pre¬ 
existing  procedures,  including  modifications  to  ensure  compatibility  with  other 
segments  of  AARC. 

The  calculation  of  nitric  oxide  radiance  observed  over  a  line  of  sight  pass¬ 
ing  through  various  regions  of  the  atmosphere,  and  possibly  through  a  zone  of 

auroral  dosing,  requires  that  AARC  divide  the  line  of  sight  into  segments.  For 
each  segment,  AARC  calculates  the  altitude,  temperature,  atmospheric  composition 
and  the  presence  and  intensity  of  auroral  NO  excitation.  These  data  allow  AARC 
to  calculate  for  each  segment  the  band  emission  rates  for  the  NO  and  N0+  funda¬ 
mental  and  overtone  bands.  These  are  then  summed  over  all  segments  to  obtain 
the  band  radiance  at  the  point  of  observation. 

AARC  then  distributes  the  band  radiances  to  all  the  rotational  lines  which 
make  up  the  band  according  to  the  emission  strength  of  the  line  as  determined 
from  the  A.F.G.L.  Atmospheric  Absorption  Line  Parameters  Compilation  <1>. 

This  set  of  line  intensities  and  positions  is  then  convolved  with  a  spectro¬ 
meter  response  function  to  generate  the  synthetic  spectrum,  which  is  output  to  a 

file  in  a  form  readable  by  a  subsequent  plotting  routine. 

AARC  allows  the  user  to  specify  all  of  the  parameters  which  characterize 
the  observation  situation.  These  include  the  position  of  the  observer  and  the 
geometry  of  the  line  of  sight,  the  character  and  boundaries  of  the  auroral  dosing, 
and  the  wavelength  range  of  the  synthetic  spectrum  along  with  the  nature  of  the 
spectrometer  response  function.  The  parameter  input  procedure  has  been  designed 
to  be  highly  user  friendly  by  presenting  the  user  a  series  of  menues ,  allowing 


the  modification  of  any  or  all  parameters.  All  parameters  default  to  the  values 
used  in  the  previous  run  of  AARC.  This  is  accomplished  by  storing  the  values  of 
the  parameters  in  a  data  file  which  is  updated  with  each  run  of  AARC. 


In  addition,  AARC  may  be  run  in  a  non-interactive  mode  which  does  not  poll 
the  user  for  input.  In  this  mode,  the  user  may  specify  input  parameters  by 
making  modifications  to  the  data  file  holding  the  defaults  for  the  input  param¬ 
eters  . 

A  complete  report  describing  the  first  release  of  AARC,  and  showing  exactly 
what  forms  are  required  for  the  input  data,  is  being  prepared  with  the  help  of 
ARCON  personnel,  and  will  be  released  early  in  1985.  The  following  sections 
provide  a  functional  outline  of  the  sections  of  the  code  which  have  been  developed 
or  substantially  modified  by  ARCON  personnel. 


2.  CALCULATION  OF  ELECTRON-ION  PAIR  PRODUCTION  PROFILE 

The  calculation  of  the  volume  emission  rate  along  the  line  of  sight  requires 
AARC  to  first  calculate  the  rate  of  electron  energy  deposition  as  a  function  of 
altitude  for  the  electron  energy  distribution  function  specified  by  the  user. 

The  energy  distribution  function  of  the  electron  flux  at  the  top  of  the 
atmosphere  may  be  specified  in  terms  of  a  Maxwellian  distribution. 


f (E)  ~  A-E*exp(-E/a)  ,  (1) 

where  f (E)  is  the  flux  of  electrons  of  energy  E  at  the  atmospheric  ceiling  in 
2 

electrons/cm  -sec-keV,  cx  is  the  characteristic  energy  in  keV,  and  A  is  a  normal¬ 
ization  constant  depending  upon  a  and  Q  ,  the  total  energy  incident  on  the  atmos- 
2  m 

phere  in  ergs/cm  -sec: 


A 


_ _ 

2-1.602  x  10_9-a3 


(2) 


a. 


To  define  this  dist 
and  the  total  energy 


ribution , 
flux,  Q^. 


the  user  must  specify  the  characteristic  energy, 
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The  electron  flux  may  also  be  specified  via  a  Gaussian  distribution. 


f (E)  =  A-E-exp 


(E-E  )2* 

_  O 

w 


(2) 


where  E  the  energy  at  the  distribution  maximum,  w  is  the  characteristic  width 
o 

of  the  distribution,  and  A  is  a  normalization  constant  determined  from  the  total 
energy  flux.  Here  the  user  must  specify  the  two  constants,  and  w,  and  also 
the  total  energy  flux. 

The  choice  of  Maxwellian  distribution  is  appropriate  for  continuous  aurora, 
while  Gaussian  distribution  is  more  appropriate  for  discrete  aurora. 


The  calculation  of  the  corresponding  energy  dosing  profile  employs  a  semi- 
empirical  method  due  to  Rees  <2>,  which  utilizes  an  experimentally  determined 
energy  dissipation  function  <3>  shown  in  Figure  33.  Here  is  the  fractional 

energy  dissipation,  z  is  the  air  mass  penetrated  by  the  monoenergetic  electron 
beam  in  gm/cm2,  and  R  is  the  maximum  range  of  the  beam  given  by  <3^: 

R  «  4.57  x  10“6  E1’75  ,  (3) 


where  E  is  the  energy  of  the  beam  in  keV.  The  energy  dissipation  function, 
is  relatively  insensitive  to  the  beam  energy  for  energies  ranging  from  0.4  keV  to 
200  keV. 


The  ion-pair  production  rate  in  pairs /cm 3-sec  is  given  by 


q(z) 


P(z) 

AE. 

ion 


f  (E) 


(4) 


where  p(z)  is  the  atmospheric  density  in  gm/cm3  at  the  altitude  z,  and  AE.  is 

ion 

the  energy  for  production  of  one  electron-ion  pair  (0.035  keV/pair). 

AARC  employs  Gaussian  quadrature  to  integrate  eq.  (4)  over  the  energy  range  of 
interest . 

Alternatively,  the  user  may  choose  from  one  of  three  program-defined  dosing 

*4" 

profiles  appropriate  for  auroras  of  IBC  classes  II,  III  and  III  .  AARC  will  then 
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Figure  33.  Normalized  dissipation  function 


use  dosing  profiles  stored  within  the  program  which  were  previously  calculated 
using  a  Maxwellian  distribution  and  parameters  given  in  Table  8.  Figure  34 
shows  the  electron-ion  pair  production  profile  for  these  three  distributions. 

3 .  CALCULATION  OF  VOLUME  EMISSION  RATE 

The  calculation  of  the  auroral  NO  infrared  volume  emission  rate  requires 

2 

the  solution  of  a  set  of  balance  equations  involving  the  production  of  N(""D) 

A  -I" 

and  N(  S)  via  collision  of  energetic  electrons  with  N  and  NO  and  the  reaction 

^24 

of  N^+  with  0^  and  0,  the  production  of  N0(v)  by  reaction  of  N(  D)  and  N(  S) 

with  0  and  emission  and  absorption  from  N0(v4l)  and  NO(v-l),  and  depletion  of 
2  2  4 

N0(v)  by  collision  with  N(  D)  and  N(  S)  and  also  by  emission  or  absorption  of 
infrared  radiation. 

AARC  does  not  explicitly  solve  these  balance  equations,  but  instead  uses  a 
set  of  production  efficiencies  which  have  been  previously  calculated  by  assuming 
an  electron-ion  pair  production  rate,  solving  the  balance  equations  for  the 
corresponding  volume  emission  rate,  and  then  dividing  the  emission  rate  by  pair 
production  rate  assumed.  These  efficiencies  are  used  as  input  to  AARC,  where  they 
are  used  to  calculate  band  emission  rates  for  an  arbitrary  dosing  rate. 

The  calculation  of  the  excitation  of  N0(v=l)  from  the  ground  state  via  earth- 
shine  in  the  quiescent  atmosphere  is  more  explicit. 

We  assume  that  the  earthshine  excitation  rate  is  constant  over  the  line-of- 
siRht  and  is  given  by: 

R,,  =  1.06  x  10  ^  excitations  N0(v=l ) /sec -mol ecu!  es  N0(v=0)  .  (5) 

E 

Relaxation  of  N0(v=l)  proceeds  via  atomic  oxygen  collision 

N0(v=l )  +  0  N0(v=0)  4  0  (Rl) 

k,  =  6.5  x  10'H 


Figure  34.  ELECTRON-ION  PAIR  PRODUCTION  RATE  (pairs /cm 3-sec) 


via  collision  with  molecular  oxygen: 


NO(v=l)  +  02  «-*-  N0(v=0)  +  02 

k2  =  3.0  x  10” 


and  also  via  spontaneous  emission. 


Under  equilibrium  conditions. 


[N0(v=l )  ] 


Re  +  k_j 10J  +  k_2[02] 
A10  +  klf^"  + 


(R2) 


(6) 


where  is  the  Einstein  coefficient  for  emission  from  N0(v=l).  The  volume 
emission  rate  is,  of  course,  just  this  times  A^q. 


it.  CALCULATION  OF  RADIANCE  AT  THE  OBSERVATION  POINT 

To  calculate  the  radiance,  AARC  must  sum  the  emission  rate,  which  depends 
upon  tue  composition  of  the  atmosphere,  temperature,  pressure,  and  auroral  energy 
deposition  rate,  over  the  entire  line  of  sight.  The  line  of  sight  passes  through 
a  range  of  altitudes  and  regions  characterized  by  presence  or  absence  of  auroral 
dosing.  Thus  AARC  must  partition  the  line  of  sight  into  segments,  and  calculate 
the  altitude  for  each  segment  so  that  the  appropriate  atmospheric  parameters  may 
be  determined  via  linear  interpolation,  and  also  calculate  the  geographic  location 
of  the  segment,  via  simple  spherical  geometry,  so  that  the  presence  of  auroral 
dosing  may  be  determined.  AARC  compares  the  coordinates  of  the  segment  with  the 
boundaries  of  the  auroral  zone,  which  the  user  may  specify  using  either  geo¬ 
graphic  or  geomagnetic  coordinates. 

The  summation  of  the  NO  infrared  emissions  begins  at  the  far  end  of  the  line 
of  sight,  where  it  passes  through  the  maximum  altitude  of  interest  specified  by 
the  user,  continues  through  the  tangent  point  to  the  near  point,  where  it  passes 
through  the  maximum  altitude  again.  The  boundaries  of  each  segment  are  chosen  so 
that  ea  ;h  segment  subtends  the  same  angle  at  earth  center.  Currently,  AARC  is 


i. 


1BC 

CLASS 

a 

(kev) 

iU 

(ergs/cm^-sec) 

II 

2.9 

12.9 

III 

5.0 

100.0 

III+ 

10.0 

400.0 

Table  8  Pair  production  profile  parameters  included 
in  AARC  for  three  classes  of  auroras. 


The  summation  over  the  line  of  sight  produces  a  set  of  band  radiances  for 
the  NO  and  N0+  fundamental  and  overtone  bands  up  to  the  vibrational  state  v'  =  15. 
In  addition,  AARC  calculates,  via  a  weighted  average  over  the  line  of  sight,  an 
effective  band  temperature  associated  with  each  band. 

The  next  step  is  to  distribute  the  band  radiances  over  all  ro-vibrational 
transitions  in  the  bands. 

The  intensity  of  a  particular  line  is  proportional  to  the  probability  that 
the  upper  state  is  occupied,  as  calculated  from  the  Boltzmann  factor  using  the 
effective  band  temperature  and  the  degeneracy  of  the  state,  and  porportional  to 
the  Einstein  spontaneous  emission  coefficient  as  calculated  from  the  AFGL  data¬ 
base  <1>.  The  calculation  of  individual  line  intensities  is  subject  to  the 
constraint  that  the  sum  of  all  lines  in  a  band  must  be  equal  to  the  total  band 
radiance. 

This  calculation  results  in  a  set  of  line  positions  and  intensities,  which 
is  then  convolved  with  a  triangular  spectrometer  function  to  generate  a  synthetic 
spectrum  over  a  range  of  wavenumbers  specified  by  the  user. 
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V.  APPENDIX 


The  source  code  for  function  ZVGT  is  given  in  the  following  pages. 


FUNCTION  Z VGT ( XX , Y ) 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


c 

c 

c 


VOIGT  PROFILE  ALGORITHM  DESIGNED  FOR  FAST  EXECUTION  IN  PROGRAM  NLTE. 
DOPPLER  PROFILE  IS  SUBSTITUTED  IN  SOME  REGIONS.  THE  METHOD  EMPLOYED 
IN  REGION  1  WAS  ORIGINALLY  OUTLINED  BY  PIERLUISSI  ET  AL, 

JQSRT  18,555(1977).  THE  ROUTINE  HAS  BEEN  TESTED  FOR  0<XX<30  AND  Y=0, 
. 000001<Y<10 .  IT  GIVES  AN  ACCURACY  OF  BETTER  THAN  .1%  FOR  XX<1.5  AND 
BETTER  THAN  .5%  EVERYWHERE  ELSE  IN  THIS  REGION.  THAT  IS,  THE 
ABSOLUTE  ERROR  IS  LESS  THAN  .001  OR  .005  TIMES  THE  CORRECT  FUNCTION 
VALUE.  FOR  Y<. 000001,  WHERE  THE  DOPPLER  PROFILE  IS  SUBSTITUTED, 

THE  FRACTIONAL  ERROR  MAY  BE  GREATER  FOR  XX>4  BUT  THE  ABSOLUTE  ERROR 
IS  VERY  SMALL. 


Y  =  A  =  SQRT(LN2 ) *ALFL/ALFD 

XX  =  ( NU-MU0 ) /ALFD  =  X/SQRT(LN2) 

AN  =  1/(N| *(2N+1) ) ,  N  =  0 , . . . , 56 

W6  =  WEIGHT  FOR  2-PT  GAUSS-HERMITE  QUADRATURE  TIMES  2/PI 

W1-W3  =  WEIGHTS  FOR  6-PT  GAUSS-HERMITE  QUADRATURE  TIMES  2/PI 

U6  =  SQUARE  OF  ZERO  FIXED  BY  2-PT  GAUSS-HERMITE  QUADRATURE 

U1-U3  =  SQUARES  OF  ZEROES  FIXED  BY  6-PT  GAUSS-HERMITE  QUADRATURE 


PISQ  =  2/SQRT(PI ) ; 

DIMENSION  AN (57) 

DATA  ( AN(N) ,N*1 , 57 )/l . 

1  . 23809523805 52E-01 , 

2  . 1068376068376E-03 , 

3  . 145038522231 5E-06 , 

4  . 8350702795147E-10 , 

5  . 2466827010264E-13 , 

6  . 4221407288807E-17 , 

7  . 45518467 589 28E-21 , 

8  . 3289 260 349176E-25 , 

9  . 1669761793417E-29, 
A  . 6180307588223E-34 , 

B  . 171885606 280 2E-38, 

C  . 3682493515461E-43 , 

D  . 6205657919638E-48 , 

E  . 8373419683872E-53 , 

F  . 9186429502399E-58 , 

G  . 8304614505929E-63 , 

H  . 6259184116949E-68 , 

I  . 3974252722665E-73 , 


SL  =  SQRT( LN2 ) ;  SLP 


, .3333333333333, 

. 4629629629630E-02 , 
.13227 51322751E-04, 
. 131 225 3296 380E-07 , 
. 5947794013638E-11 , 
. 14483 26 4 6436 0E-: 4 , 
.2107855191442 E- 18, 
. 1977064753878E-22 , 
. 1264107898899E-26 , 
. 5754191643982E-31 , 
. 1930357 208815E-35 , 
. 4908923964523E-40 , 
.  9687280238871E-45 , 
.  1513107949 541E-49 , 
. 1902 541 227 290E-54 , 
.  1954102 582324E-59 , 
. 1660580513451E-64 , 
. 118076183891 2E-69 , 
. 7095717391818E-75 , 


=  SQRT( LN2/PI ) 


.10, 

-7575757575758E-03 , 
. 1458916900093E-05 , 
. 1089 222 10371 5E-08 , 
. 3955429516459E-12 , 
. 80327 35012416E-16 , 
.1002516493491E-19, 
•8230149299214E-24, 
. 4 67848351 55 19E-28, 
.1916942862110E-32, 
. 5846755007469E-37 , 
. 1 36304 1261779E-41 , 
. 2483069097455E-46, 
. 3601579309810E-51 , 
. 4 2257897 54194E-56 , 
.4070135277853E-61, 
. 3 255 39546201 3E-66 , 
. 21862104 2 2954E-7 1 , 
. 1 24466 5977389E-76/ 


DATA  SLP  /. 4697186393498/ 
DATA  PISQ, SL  /l . 128379167096 , 
DATA  W1,W2,W2/. 4613135279626, 
DATA  U1,U2,U3/. 1901635091935, 
DATA  W6  >.5642150484/ 

DATA  U6  / .  5/ 


.8325546111577/ 

.09999216171035, . 2883893874868E-02/ 
1.784492748543,  5.525343742263/ 


X  =  XX*SL 

IF(XX.GT. 5 . 9 )tf 3  TO  120 

IF( Y. LT. 1 . OE-06 )GO  TO  130 

IF( Y.GT. 1 . 5 )GC  TO  105 

IF(Y.LT. .0005.AND.XX.LT. 2. )GO  TO  130 


REGION  ONE 


S  =  X*X  -  Y* Y 
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I  T  -  2.0*X*Y 

SER  =  Y  y 

SEI  »  -X 

FPR  »  Y 

•  PPI  ■  -X 

C 

NMAX  *  INT(8.7*XX)  -  4  ♦  6*INT(Y) 

N  *  8  +  4*INT( Y) 

IF ( NMAX. LT.N) NMAX  =*  N 
C 

DO  100  N  >  1,NMAX 
FNR  =  FPR*S  -  FPI*T 

FNI  =  FPI*S  +  FPR*T 

SER  -  SER  +  FNR*AN{N+1) 

SEI  =  SEI  +  FNI*AN(N+1) 

FPR  =  FNR 

100  FPI  *  FNI 

ZVGT  =  SLP*EXP ( -S ) * { OS ( T ) * ( 1 . 0  -  PISQ*SER)  -  S7N(T)*PISQ*SEI) 

RETURN 

C 

C  REGION  TWO:  6-PT  GAUSS-HERMITE  QUADRATURE;  Y>1.5 
C 

105  S  -  X*X  -  Y*Y 
T  *  2.0*X*Y 
T2  *  T*T 
T  ■  T*X 
FI  *  S  -  U1 
F2  ■  S  -  U2 
F3  =*  S  -  U3 

ZVGT  *  SLP*(W1*<T  -  Y*F1)/(F1*F1  ♦  T2)  + 

1  W2*(T  -  Y*F2)/<F2*F2  +  T2)  + 

2  W3*(T  -  Y*F3)/(F3*F3  +  T2)) 

C 

RETURN 

C 

C  REGION  THREE:  2-PT  GAUSS-HERMITE  QUADRATURE  (FATER  THAN  DOPP);  XX>5. 
C 

120  T  =  2 . 0*X*Y 

SEI  =  T*T  1 

T  *  T*X 

FPR  *  X*X  -  Y*Y  -  U6 

ZVGT  *  SLP*W6*(T  -  Y*FPR)/(FPR*FPR  +  SEI) 

RETURN 

C 

C  REGION  FOUR:  SUBSTITUTE  DOPPLER  PROFILE 
C 

130  ZVGT  «  SLP*EXP(-X*X) 

RETURN 

END 
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